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Preface 


The  main  problem  scheduled  for  the  year  19^3  by  the  team  v^'orking 
on  VLP  problems  since  December  196O  is  the  explanation  and 
elimination  of  discrepancies  between  the  results  of  theoretical 
calculations  and  experimental  measurements  that  still  existed. 

For  this  purpose  all  factors  influencing  the  propagation  of 
VLP  waves  had  to  be  studied  theoretically  and  the  experimental 
methods  of  measurement  had  to  be  further  developed.  The  improvement 
of  all  apparatus  used  so  far  as  well  as  the  expansion  of  measure¬ 
ments  over  much  greater  distances  at  v/hich  clear  effects  are  to 
be  expected  according  to  theoretical  considerations,  are  of  special 
importance.  Measurements  v/ere  conducted  in  German  and  Dutch  mines 
to  explain  the  effect  of  different  geological  bodies  on  wave 
propagation}  furthermore,  a  number  of  conductivity  measurements 
were  made.  V/e  mainly  thank  our  "landlord'*  Bergrat  Dipl.  Ing.  Frottner, 
director  of  the  Montan  "'erke  Brixlegg  for  making  it  possible  for 
us  to  work  and  move  about  freely  with  our  measuring  equipment  in 
the  mines  of  Schwaz  and  St.  Gertraudi. 

He  also  gave  us  permission  to  install  our  mine  lab,  and  some  of 
his  miners  did  the  necessary  blasting.  again  want  to  express 
our  best  thanks. 

The  directors  of  the  Dutch  coal  pit  Oranje  Nassau  III  m.-ule  it 
possible  for  us  (by  special  safety  measures  and  by  unsupported 
galleries)  to  conduct  an  interesting  measurement  along  a  coal  seam. 

Fe  want  to  thank  them  very  much  for  their  great  efforts. 

Furthermore  I  want  to  express  my  special  thanks  to  the  directors  of 
the  Salzgitter  ore  mine,  Vereinigte  Kaliv/erke  Saldetfurth,  Siegerland 
ore  mine  and  all  those  who  assisted  us  in  our  subsurface  measurements. 


Innsbruck,  '^-ecember  19^3 


V/olfram  Bitterlich 
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of  a  nev;  VLP  receiver,  A  completely  new  way  was  chosen  to  build 
antenna  circuits,  amplifier  circuits,  and  selective  circuits 
from  different  vievpoints  and  by  means  of  nuvistors,  (See  Fig.  III). 
The  first  results  which  are  promising  show  a  progress  in  our 
efforts  to  build  an  optimum  receiving  device.  The  sensitivity, 
band  width,  and  the  signal-to-noise  ratio  were  improved. 

These  studies  will  probably  be  published  in  I964,  as  the 
development  is  still  being  continued. 

An  easily  portable  measuring  device  v;as  built  and  numerous 
experiments  in  different  geological  bodies  were  made  in  order 
to  determine  the  electrical  data  of  the  rock,  especially  the 
conductivity. 

In  '•963,  45  expeditions  to  mines  and  underground  measurements 
were  made, 

1.2  Field  strength  measuring  device 

The  -field  strength  measuring  device  built  in  19^2  had  to  be  improved 
in  various  respects.  In  principle,  however,  the  construction 
remained  the  same.  The  following  items  have  to  be  pointed  out: 

(1.)  Redesign  of  the  Input  stage;  lov/ering  of  the  noise  level,  and 

at  the  same  time  increase  of  amplification,  A  limiter  circuit  at 
the  input  prevents  this  stage  to  be  overdriven. 

(2.)  The  selectivity  was  increased  by  building  in  a  helical 

potentioiiiGter ,  simultaneously,  the  amplification  was  made  more 
constant  (better  than  +  2  db  v^hen  txining  over  a  frequency  range 
of  1  :  10).  The  frequency  range  of  50  cps  -  3P  kc/sec  can  now 
be  covered  continuously  in  three  ranges. 

(3.)  The  use  of  a  millivol tmeter  proved  to  be  very  suitable  so  that 
tlie  antenna  voltage  on  the  amplifier  input  in  this  developed 
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device  can  be  read  directly  from  the  instrument.  The  scale  is 
calibrated  in  microvolts  or  millivolts.  Voltages  bet7?een  0.1  p' 
and  100  mv  can  be  measured. 

(4.)  Special  attention  ^vas  paid  to  the  temperature  stability  of  the 
device,  as  the  environmental  temperature  varied  from  4°C  to 
35°C.  The  circuit  diagrams  of  Figs.  1,1  to  1,4  show  that 
considerable  independence  of  temperature  and  fluctuations  of 
the  operating  voltage  was  reached  by  means  of  directly  coupled 
amplifier  stages  and  sufficiently  dimensioned  negative  feed¬ 
backs.  The  opierating  voltage  is  kept  constant  by  an  electronic 
circuit,  battery  voltage  fluctuations  from  20  v  to  14  v  having 
no  effect  on  the  measuring  accuracy. 

(5.)  A  special  calibrated  oscillator  t^s  provided  to  make  it 
possible  that  the  device  can  be  recalibrated  any  time.  It 
produces  rectangular  voltayc  pulses  of  approximately  1,3  kc/seo 
whose  value  is  made  independent  of  the  environmental  temperature 
by  means  of  careful  stabilization. 

This  calibrated  voltaye  is  fed  to  the  amplifier  input  in 
"calibration"  position.  The  indicator  of  the  instrument  can 
be  adjusted  to  the  corresponaing  calibration  value  by  means 
of  the  potentiometer  . 

(61)  In  principle,  the  superheterodyne  unit  shov'n  in  Fig,  1,4 

remained  unchanged.  The  only  dilference  occurs  in  the  ou\put 
amplifier  which  nor' contains  complementary  transistors  in  order 
to  avoid  the  use  of  any  transformers  v/hose  magnetic  stray  field 
might  affect  the  antenna  stages. 

The  field  strength  meter  hunt  according  to  those  principles 
is  shown  in  Fig.  IV.  The  gr 'at  advantage  of  this  dev ' ce  is 
its  reliable  constru.ction  and  its  low  weight. 
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1,5  The  amplifiers  of  the  transmitter  and  the  trans..iitting 
antennas 

In  principle,  the  battery-operated  device  rpinai;.«d  unchanged, 

but  an  autotransformer  made  it  ^ossible  to  reach  a  higher 

output  pov/er.  Fig.  1,5  sho^vs  that  the  antenna  circuit  can 

thus  be  directly  connected  v/ith  the  output  transistors, 

v.’ithout  involving  any  additivial  drawbacks.  The  operating 

reliability  was  improved  by  using  an  overvoltage  fuse  for 

output  transistors.  This  is  to  avoid  that  short-term  voltage 

peaks  ■'.’rhich  might  occur  e.g,  by  tuning  to  resonance  by  means 

of  the  capacitor  dcoa Ic,  destroy  the  power  transistors. 

This  protective  circuit  consists  cf  heavy-duty  diodes  and  a 

sufficiently  dimensioned  buffer  capacitor  (11  raF). 

As  the  capacitor  decade  in  the  improved  amplifier  of  the 

transmitter  had  been  made  stronger  so  that  resonance  voltages 

up  to  2000  V  could  occur  without  any  damage,  the  mechanical 

construction  of  the  Jevice  had  also  to  be  changed 

correspondingly.  The  device  was  used  mainly  for  measurements 

in  Germany  and  Holland.  An  amplifier  of  1000  w  was  installed 

iti  the  St.  Gertraudi  mine  in  order  to  transmit  over  larger 

distances.  The  device  delivered  by  the  firm  Bryan  SAVAGE  is 

mains  supplied.  The  frequency  range  of  the  amplifier  at  full 

output  pov/er  reaches  20  ops  -  10  kc/sec.  Especially  in 

connection  with  the  transmitter  antenna  SA  VIII  (which  will 

be  described  below)  with  a  winding  surface  of  approximately 
2 

480  m  ,  distance  of  more  than  5000  n  could  be  bridged  v/ithout 
any  difficulties- 

As  wG  have  mentioned  before,  this  transmitter  was  set  up 
in  r  special  mine  laboratory  This  was  necessary  because  of 
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the  available  power  and  also  because  such  a  device  has  to  be 
set  up  under  certain  climatic  conditions.  The  atmospheric 
humidity  had  to  be  reduced  from  ^QOfo  to  standard  values,  and 
the  installations  had  to  be  made  v/ith  special  care  for  such 
pov’ers.  Fig.  I  shows  the  transmitter  and  the  mine  lab.  Two 
separate  antennas  were  fed  with  this  device. 

1.  Pivoted  frame;  SA  VII. 

In  order  to  make  possible  measurements  under  different  angles 

the  existing  frame  in  the  *'Bunte  Kluft",  with  a  winding  surface 
2 

of  500  m  ,  had  to  be  reconstructed  so  that  it  can  be  used  for 
the  provided  power  of  1  kw.  For  this  purpose,  it  was  necessary 
to  insulate  the  individual  turns  of  the  frame  sufficiently. 
Voltages  that  occur  between  the  individual  turns  or  also  with 
respect  to  ground  exceed  2000  v,  so  that  the  plastics  insulation 
of  the  wires  becomes  insufficient.  Therefore,  the  individual 
’•Tires  were  installed  at  a  certain  distance  in  oi'der  to  avoid 
a  direct  contact  of  the  wires.  Fig.  II  shows  the  new  antenna 
whose  dimensions  on  the  "'hole  remained  the  same,  namely 
5.8  m  X  4  m  coil  area. 

2.  Fixed  frame;  SA  VIII. 

For  measurements  over  larger  distances,  a  pivoted  antenna 

(SA  VII,  see  above)  coulJ  not  be  used.  The  next  step  was 

therefore  the  construction  of  a  threefold  frame  vdth  a  "inding 
2 

area  of  480  ra  .  The  dimensions  of  this  frame  are  approximately 
4'J  r  X  40  m.  The  antenna  was  laid  through  vertical  shafts  and 
the  connecting  galleries  above  en  i  belov/.  In  this  case  it  was 
also  important  to  insula tf>  '^ntenna  wire  perfectly  from  the 
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ground,  otherwise  a  considerable  portion  of  the  output  power 
would  have  been  lost. 

/  2  \ 

Sufficient  dimensioning  of  the  wire  cross  section  (2*5  ^  Gu; 
fillovtd  ug  to  choo8«  the  r«8lgt*ncg  of  tha  so  saalX 

that  a  current  of  more  than  25  a  passes  throu^  the  antenna  at 
full  output  power.  The  two  transmitter  antennas  differ  by  a 
factor  which  can  be  expressed  by  winding  area  times  antenna 
current.  The  following  values  are  obtained: 

SA  VII  :  750  a*m^ 

SA  VIII  :  1350  a-m^ 

SA  VIII  is  thus  superior  to  SA  VII  by  the  factor  1,8.  The  only 
drawback  is  the  fact  that  SA  VIII  is  not  pivoted  so  that  it 
can  only  be  used  for  measurements  at  one  single  angle  ■5'. 

,4  Mine  laboratory 

A  special  room  with  a  power  supply  (580  v  50  cps)  had  to  be 
equipped  for  setting  up  the  large  mains-supplied  transmitter. 

As  it  was  desirable  that  a  number  of  measuring  instruments  should 
be  left  permanently  in  the  new  lab,  the  atmospheric  humidity  in 
the  mine  had  to  be  reduced  from  ^OOfo  to  standard  values. 

This  is  done  by  tvfo  powerful  air  drying  apparatus .  In  order  to 
increase  their  efficiency,  the  room  is  permanently  kept  at  a 
temperature  of  more  than  25®C,  It  is  very  difficult  to  air- 
condition  the  laboratory  because  the  surrounding  rock  is  very 
damp  and  the  floor  of  the  gallery  is  covered  \7itn  percolating 
water. 

The  mine  lab  shov/n  in  Fig.  I  is  equipped  v;ith  the  following 


devices : 
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(1)  Low-frequency  amplifier  ranging  from  20  cps  to  10  kc/sec 
with  1  kw  output  power. 

(2)  High  frequency  amplifier  for  frequencies  ranging  from 
100  kc/sec  to  approximately  600  kc/sec  v^ith  a  maximum 
output  power  of  1200  w. 

This  transmitter  is  intended  to  be  used  for  conductivity 
measurements.  The  theoretical  bases  of  these  measure¬ 
ments  are  outlined  in  our  second  annual  report  [1]  and 
in  [4]. 

(3)  A  number  of  instruments  such  as  an  oscilloscope,  a 
multirange  voltmeter,  an  RC  generator,  an  AP  amplifier, 
a  VLP  receiver,  and  a  tube  voltmeter. 

The  laboratory  will  become  very  important  for  future  work, 
because  measurements  of  VLP  transmitters  and  atmospherics  over 
extended  periods  of  time  are  planned.  The  position  of  the  lab 
in  the  mine  is  shown  in  Fig.  1,21  (a).  This  site  was  chosen 
because  it  offers  the  possibility  of  transporting  heavy  instru¬ 
ments  with  the  mine  railv/ay  right  to  the  mine  lab,  and  because 
antennas  can  easily  be  set  up  in  its  immediate  neighborhood. 
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Propagation  measurements  of  VLF  waves 

As  \7e  have  briefly  mentioned  in  the  preface,  the  experimental 
studies  of  19^3  concentrated  on  two  points: 

1 .  The  most  important  tasks  were  studies  of  propagation 
in  geological  media  which  differ  largely  from  those  examined 
so  far.  The  excursion  to  Northern  Germany  and  Holland  yielded 
a  number  of  interesting  results  concerning  this  sphere  of 
problems.  We  conducted  measurements  in  five  different  mines 

in  media  of  high  conductivity  and  also  in  well  insulating  rock. 

The  possibility  of  measuring  under  simple  geometrical  conditions, 
i.e.  conditions  under  which  the  earth's  surface  can  actually  be 
looked  upon  as  a  plane,  was  of  great  interest.  Our  theoretical 
calculations  are  all  based  on  this  simplified  assumption  so  that 
the  measurements  made  in  the  Tyrol  are  only  approximately 
comparable  to  the  results  obtained  theoretically. 

2.  Measurements  over  very  large  distances,  i.e.  more  than 
one  kilometer.  They  were  expected  to  yield  results  concerning 
our  apparatus  and  especially  the  construction  of  devices  that 
would  make  possible  transmission  over  still  greater  distances. 

At  the  same  time  another  problem  could  be  solved,  namely  the  problem 
whether  the  sky  wave  or  the  ground  wave  is  responsible  for  the 
propagation.  This  problem  is  to  be  discussed  in  part 

1,,5  Salzgitter  mine  Konrad  I  (Germany) 

The  measurements  conducted  in  the  iron  ore  mine  Konrad  I  yielded  e 
number  of  very  interesting  results,  as  the  conductivity  of  the 
studied  geological  medium  differs  considerably  from  those  studied 


so  far. 
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The  experimental  results  were  comfeired  v/ith  the  theory  given 
in  fl]  that  holds  for  a  dipole  in  an  unbounded  homogeneous 
medium.  As  the  conductivity  vas  unexpectedly  high,  the  numerical 
data  given  in  {ij  had  to  be  calculated  again  for  the  special 
values  of  conductivity  and  distance. 


1 .  Theoretical  considerations 


According  tt  'e  theory  developed  in  the  second  Annual  VLP  Report 
(1962),  which  holds  for  a  magnetic  dipole  in  an  unbounded  and 
conducting  medium,  the  mag:ietic  field  strength  was  calculated 
as  dependent  on  the  distance  r,  the  angle  and  the 
conductivity  0.  Below,  the  quantity  (see  p.l  1  expression  40 
and  41  in  the  above  report)  instead  of  the  conductivity  a 
shall  be  considered  to  make  the  calculation  more  general, 
otherwise  the  results  "'ould  be  bound  to  a  certain  measuring 
frequency,  but  this  is  to  be  avoided.  For  this  purpose,  the 
function  h(r,'8‘,A)  (expression  52  of  the  report)  was  calcula,ted 
for  a  number  of  distances  and  conductivities. 


The  function 


h(r,,'>,A)  =  (  1  ^  24  ^  2,^2 


2  '  2  r 

r  ^ 


.  2  a  /  1  2/5?  4^  .  2  X 

+  sin  (-7  +  -v  —5-  +  4^  ) 


)  + 


(1) 


4  3  2  r 

r  r  r 
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is  calculated  in  Table  I  for  the  values  A  =  4 •86*10  ; 

1.085*10"^;  3. 45*10''^;  1.09*10“\  For  the 

frequency  5  kc/sec  used  for  ese  "‘•asurements ,  the 

-2  -1  -1 

conductivity  would  le  r  ^  10  •  10  ;  2*  10  ; 
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1  ‘t 

l(i}.  ~  m~  ) .  The  graphical  representations  shown  in  Figs.  1,6 
through  1 ,9  give  a  good  illustration  of  the  directional  pattern 
as  dependent  on  the  distance  at  constant  conductivity  on  the  one 
hand,  and  at  a  given  distance  as  dependent  on  the  conductivity 
on  the  other  hand.  For  a  better  survey,  the  value  of  the 
function  h(r,'5',A)  was  plotted  in  the  diagi’am  in  radial  direction 
(logarithmic  scale).  Fig.  1,8  clearly  shows  the  result 
obtained  previously,  namely  that  the  field  strength  component 
H  (in  the  direction t  =  90°)  ^  certain  conductivity  has  a 

maximum  at  a  given  distance  and  frequency. 


The  definition  of  the  function  (p(r,/l)  in  the  form 

,0 


/  !Hl  ^  =  0 

- 

iHl  900 


(2) 


proved  advantageous  for  further  studies.  Table  II  gives  some 
computed  values  of  this  function  v/hich  are  plotted  in  Fig,  1,10. 
By  means  of  this  function  which  on  the  whole  is  consistent  \7itVt 
the  function  V(r)  in  the  report  (see  p.  l6,  expression  55)  >  l^he 
conductivity  can  be  determined  from  the  value  of  (p(r,A)  v/hich 
can  easily  be  measured  by  experiment,  and  from  Pig,  1,10  via 
The  considerable  decrease  of  the  value  of  the  function 
characterizing  the  far  field  of  the  dipole  indicates  the 
predominance  of  the  component  H^. 

In  the  present  report  it  shall  not  be  discussed  hov;  far  the 
boundary  of  the  medium  can  be  neglected  and  how  great  the  error 
thus  becomes,  since  the  existing  results  at  present  are  being 
compared  with  the  experimental  data. 

2.  Measurements  in  the  Salzgitter  mine  (Germany)-  shaft  region 
Konrad  I 

The  geological  structure  of  this  region  is  characterized  by  a 


{ lx, if  A) 
l*Jkd 
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variety  of  flat  strata  formed  by  repeated  flooding  by  the  sea.  The 
exploitable  iron  ore  field  of  the  Upper  Jurassic  is  accompanied  by 
zones  of  marly  limestone.  The  rock  on  top  of  the  mine  has  an 
average  height  of  1000  m,  the  inclination  of  the  deposit  to  the 
horizon  being  22°.  The  concentration  of  this  oolitic  ore  is 
approximately  JO^,  the  average  size  of  the  individual  brown 
iron  grains  is  0*5  mm. 

Suitable  measuring  sites  with  specific  conditions  in  the  mine  itself 
were  not  easy  to  find.  Regions  with  v/ires  of  the  mining  railway, 
power  cables,  or  even  iron  formv/ork  had  to  be  avoided.  Two  short 
drifts  18  and  25  ni  long  in  the  1000m  bed  seemed  to  be  most  suitable. 
Despite  the  short  distance  betv/cen  transmitter  and  receiver, 
characteristic  antenna  patterns  could  be  measured  owing  to  the 
high  conductivity  of  the  ore  body  -  a  d-c  measurement  yielded 
o  =  2*10'’  (iTl  m  ).  The  working  method  for  finding  these  diagrams 
had  already  been  described  in  detail  in  the  above  annual  report. 

The  measurements  with  the  receiving  antenna  were  conducted  such 
that  the  antenna  is  revolved  3o0°  about  its  vertical  axis  at  a 
constant  angle  v'.  This  rotation  corresponds  to  a  change  in  angle 
between  the  dipole  axis  of  the  receiving  antenna  and  the  field 

vl 

direction.  The  angle  v"  is  the  angle  betv/cen  the  direction  of  the 

dipole  axis  of  the  transmitter  antenna  and  the  direction  to  the 

receiving  point.  The  measured  vali’ns  obtained  by  this  method  are 

given  under  U^  in  Tables  III  and  IV,  Column  v'as  obtained  from 

these  values  by  averaging  the  measured  values  belonging  to  4^,  4>+180° 

4; +90°.  The  mean  values  thus  obtained  are  shown  in  Pigs.  1,11  and 

1,12.  The  third  column  gives  the  mean  values  calculated  bv  the 

method  of  least  squares,  v/ith  allowance  for  the  dependence 

(j  =  Uo  ‘Cosij;.  According  to  Pigs.  1,11  and  1,12,  this  dependence 
e  '-'O 

on  the  angle  4^  obviously  doc:;  .'"ol  occur,  since  the  minimum  of  the 
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curve  does  not  assume  the  value  zero,  but  decreases  only  down  to  a 
certain  fraction  of  the  maximum  value  (see  above).  In  this  case, 
the  situation  is  probably  more  complex  than  it  was  shown  by  the 
above  theory. 

The  voltage  induced  in  the  antenna  calculated  according  to  the  law 
of  induction  is:  *  d^  *  ^  df.cos(J;  *  ^B»P»cos4>  (3)» 

If  B  is  considered  to  be  a  constant  vector  which  in  the  general  case 
may  be  looked  upcn  as  a  resultant  of  the  primary  field  and  a  series 
of  secondary  fields,  then  it  follows  from  the  above  equation  that  the 
induced  voltage  becomes  zero  for  the  angle  4j  *  90°.  Hence,  no  con¬ 
clusions  can  .be  made.  from;  the  measured  deviation  to  the  disturbing 
effect  of  secondary  fields  caused  by  the  boundary  of  the  ore  body. 

The  most  obvious  assumption  that 
can  still  explain  the  measurement 
is  that  of  a  highly  inhomogeneous 
field.  In  this  case,  the  situation 
at  the  receiving  antenna  can  thus  be 
represented. 

The  induced  voltage  can  be  represented  by  the  expression  (4)*  The 
angles  4^  and  (|»2  are  different  owing  to  the  field  ihhomogeneity  and 
are  interrelated  (see  diagram)  in  the,  following  ways 

®2*  °°®‘^2  +  4'2  +  «  »  ‘'80^  (5) 

This  relation  can  be  substituted  in  (4)>  ot  being  a  measure  of  the 
inhomogeneity.  Prom  these  assumption  it  follows  that  the 
induced  voltage  need  not  become  zero  for  whatever  value  of 
if  the  receiving  antenna  is  turned,  i.e.,  if  4^  is 
changed.  This  is  consistent  with  the  actual  observation. 

Here,  the  effect  is  explained  only  qualitatively,  not  con¬ 
sidering  the  reasons  for  the  great  influences  cn  the  trans- 
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nntter  field.  The  following  values  of  the  function  (p(r,^)  for 
18  and  25  ni  are  calculated  from  the  two  measurements: 


r(m) 

1.15 

18 

1 

0.75 

25 

These  values  are  plotted  in  Fig.  'tjfO.The  value  obtained  with  good 
approximation  for  the  conductivity  by  interpolation  is 
a  =  5*10”  (Q”  ’m  ).  The  error  is  no  more  than  25^  which 
considering  the  small  distance  is  very  good.  Furthermore,  this 
conductivity  is  also  consistent  with  that  measured  T;ith  direct 
current:  2 •  10'"^(52'  ^  ^ra  '), 

For  a  better  understanding  of  Figs .  it  must  be  added  that 

the  envelope  of  all  curves  gives  the  field  strength  diagrams 
represented  in  Pigs.  1,fcthrough  .  The  conducted  measurements  shovi 

that  ip(r,6)<1  for  a  distance  of  25  m,  hence  measurements  for 
distances  of  50  m  in  the  studied  rock  are  already  conducted  in 
the  far  field.  This  fact  is  of  special  interest,  since  the 
conditions  in  the  mines  of  Schwas  and  St.  Gertraudi  are  completely 
different;  here,  <p(r,6)  '^as  found  to  decrease  but  slightly  even  at 
1000  m.  For  this  reason,  detailed  studies  with  allowance  for  the 
upper  and  lower  boundaries  of  the  ore  body  can  be  expected  to 
yield  interesting  results, 

1,6  Potassium  mine  Hansa  III,  Empelde  near  Hannover  (Germany) 

In  Konrad  I,  the  conductivity  of  rock  was  generally  three  orders 
of  magnitude  higher  than  in  the  mines  examined  so  far,  in  the  mine 
Hansa  III,  however,  measure^'ntc  were  made  in  a  medium  whose 
conductivity  is  so  lov:  that  it  can  approximately  be  consid'^red 
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an  insulator.  The  measurements  were  made  in  a  large  salt  block, 
under  conditions  which,  owing  to  the  low  humidity,  suggested  the 
results  to  be  similar  to  those  measured  in  air. 

The  bottom  of  a  bed  1000  m  below  ground  was  chosen  as  meascring 
site,  because  contained  no  rails  and  no  cables;  furthermore, 
we  used  a  number  of  gangways  far  away  from  power  cables. 

Despite  the  favorable  measuring  conditions,  a  number  of  unexpected 
interferences  occurred. 

(1)  Owing  to  the  poor  damping  in  salt  bodies,  considerable 
interferences  by  the  underground  railways  in  higher  gangways 
were  observed. 

(2)  The  salt  is  being  poured  through  inclined  shafts. 

Owing  to  the  low  content  of  humidity,  it  gets  electrostatically 
charged,  Continuoun  discharges,  however,  cause  considerable  and 
irregular  fluctuations  in  deflection  of  the  measuring  instrument, 
making  exact  measurements  very  difficult. 

For  the  planned  measurements,  the  transmitter  was  set  up  at 
PP  107  (see  Pig.  1,15)  the  antenna  being  aligned  such  that  the 
angular  deflections  for  the  individual  points  of  measurement 
remained  as  small  as  possible.  Table  V  gives  the  measured 
values  for  the  individual  points.  The  different  '0  values  were 
not  taken  into  account  in  the  graphical  representation  of  Pig.  1,2-. 

The  magnetic  field  strength  is  plotted  in  arbitrary  units  so 
that  a  larger  number  of  measurements  of  several  mines  can  be 
represented  in  one  figure.  Furthermore,  the  individual  curves 
were  shifted  upward  or  downward  in  parallel;  hence,  conclusicns 
from  the  position  of  two  curves  must  not  be  made. 
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The  dependence  h(r,-d‘,(i>)  is  given  in  Table  IV  and  Fig.  1,14 
with  the  angles  d"  and  4-  being  changed  at  a  constant  distance 
of  240  m.  From  this  measurement  the  value  of  ip(r,A)  :7as 
determined  and  found  to  be  approximately  2.  This  result  also 
suggests  that  the  conductivity  of  the  medium  is  low. 

Although  damping  in  the  salt  body  was  poor,  the  transmitter  GBR 
which  is  very  strong  above  ground  could  not  be  received  anywhere 
in  the  mine.  Reception  measurements  were  conducted  200  -  1000  in 
below  ground.  There  are  probably  two  rea;^  ns  for  this  negative 
result  of  the  measurements : 

(1)  There  is  a  layer  of  ground  water  at  a  depth  of  approx¬ 
imately  180  m  which  contains  much  salt  and  thus  screens  signals 
from  above  ground  completely,  no  matter  whether  they  come  from 
VLP  transmitters  or  atmospherics. 

(2)  The  propagation  mechanism  of  VLF  signals  is  such  that 
there  exists  a  considerable  difference  whether  the  wave  penetrates 
the  earth's  surface  vertically  or  horizontally.  The  second 
viewpoint  is  supported  by  a  number  of  measurements  in  other  German 
mines  where  there  is  certainly  no  groundwater  layer  that  contains 
salt,  and  also  by  good  receiving  conditions  in  the  mines  of 
Schwaz,  Lafatsch,  and  Gertraudi,  where  the  receiving  conditions 
were  very  good  in  similar  or  greater  depths. 

1,7  Ore  mines  of  Fuesseberg  and  Georg  (Siegerland  -  Germany) 

The  conditions  in  these  two  mines  were  similar  to  those  in  the 
mines  of  Schwaz  and  Gertraudi.  Iron  ore  (siderite)  which  is 

/%T»m  a  o4-oc»yn  a  r\r\T*A  v  1  m  a  ^  o  T  V 

10  m  thick.  In  Fuesseberg  it  is  straight,  whereas  in  the  Georg 
mine  it  is  shaped  like  a  horseshoe  (see  Figs.  1,15  and  1,1?) 
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The  measuring  sites  were  chosen  in  such  a  way  that  most  cf  the 
line  connecting  transmitter  and  receiver  passed  through  ore, 
or  such  that  the  points  of  measurement  were  in  front  of  or 
behind  the  plate.  This  arrangement  made  it  possible  that  the 
effect  of  an  inhomogeneous  medium  on  the  propagation  of  VLF  waven 
could  be  studied.  The  curves  that  belong  to  Table  VII  and 
Table  VIII  are  shown  in  Fig.  1,20.  On  the  whole,  these  curves 
have  the  same  course  as  those  of  our  measurements  in  Schwas  [1] 
which  indicates  similar  conductivity  of  the  rock.  Conductivity 
measurements  w5 th  direct  current  7/hich  still  have  to  be 
discussed  show  that  the  dc?i  rock  and  the  ore  in  the  tv/o  mines 
are  of  almost  equal  resistivity.  Hence,  a  small  difference 
between  the  measurements  in  front;  of  and  behind  the  ore  plate 
is  observed.  Tables  IX  and  X  give  the  values  of  another 
measurement  made  under  the  following  conditions: 

The  transmitter  tvas  placed  at  in  a  gangway  100  m  higher, 
the  axis  of  the  transmitter  antenna  being  horizontal  so  that 
approximately  equal  angles  are  obtained  for  all  points  of 
measurement.  Then,  the  field  strength  was  measured  at  the  given 
points,  the  receiving  antenna  having  a  position  that  corresponds 
to  •  0°  (i.e.  maximum  input  voltage).  This  measurement  could 
not  be  represented  graphically,  as  the  measured  values  depend  on 
two  independent  quantities,  the  uistance  r  and  the  angle  The 
measurements  given  in  the  present  report  are  mainly  to  bo  used 
for  comparing  the  experimental  with  the  theoretical  values  of  a 
dipole  at  a  certain  depth  below  the  earth’s  surface  i?!*  Owing 
to  the  complicated  arrar  ement,  the  evaluations  have  not  been 
completed  as  yet.  In  the  discussed  c.'se,  the  distance  between 
transmitter  and  the  earth's  surface  was  approximately  equal  to 
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the  distance  between  transmitter  and  receiver  so  that  a 
noticeable  deviation  of  the  measured  value  from  the  theoretical 
value  of  a  dipole  in  an  unbounded  medium  was  to  be  expected 
according  to  f  7 ] • 

As  in  the  case  of  Konrad  I,  the  directivity  patterns  could  be 
plotted  for  several  distances  for  both  mines.  Different  from 
the  measurements  in  Konrad  I,  the  curves  here  show  a  clear 
minimum  (see  Tables  XI  -  ’'III,  column  U  ).  The  effect  of  the 
bounds  of  the  ore  plate  is  very  small,  because  the  electrical 
properties  of  the  two  media  differed  but  slightly.  The 
graphical  representation  of  the  measurements  is  given  in 
Figs.  1,16;  1,18;  1,19  balanced  values  being  plotted 

radially  in  linear  scale. 

1,6  Oranje  Nassau  mine  in  Heerlen.  Tolland 

"Nassau  III"  offered  us  the  only  possibility  of  measuring  in 
a  coal  mine.  Owing  to  the  way  the  coal  is  deposited  there,  we 
had  to  measure  in  the  direction  of  a  seam  approximately  2  m 
thick.  The  management  of  the  mine  had  ordered  that  timber- 
propped  galleries  were  driven  into  the  rock  so  that  there  was 
no  iron  in  the  i^'mediate  neighborhood  of  transmitter  and 
receiver.  The  two  receiving  sites  were  placed  at  a  distance  of 
500  and  500  m  from  the  transmitter.  The  transmitter  antenna 
had  an  angle  of  •>  =  0°  and  '' =  90”  with  respect  to  the  receiving 
point,  '"e  were  surprised  to  find  that  our  conductivity  measure¬ 
ments  vnth  direct  ciirrent  showed  the  coal  to  be  of  poor  conduc¬ 
tivity,  whereas  the  rock  arove  and  below  it  is  highly  conductive. 
The  next  chapter  will  give  the  values  of  d-c  measurements  which 


show  a  great  difference  bet/^^n  -roal  and  dead  rock. 
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propagation  in  this  case  took  place  along  a  thin  layer  of  poor 
conductivity,  bounded  by  tv/o  media  of  high  conductivity. 

A  geometrical  arrangement  of  this  type  which  in  principle  is 
found  in  everv  coal  mine,  seems  to  be  suitable  for  the 
propagation  of  VLF  waves  according  to  the  results  of  the 
measurements.  We  have  therefore  planned  theoretical  calcula¬ 
tions  for  studying  the  propagation  in  a  "sandwich". 

Unfortunately,  these  probleris  cculd  net  be  sufficiently 
studied  on  the  spot,  but  ihe  fcl] owing  factors  seem  to  be  of 
special  importance  for  the  propagation  along  a  coal  seam: 

(1)  Thickness  of  the  seam  ?s  compared  to  the  total 
distance  between  the  transmitter  and  receiver.  The  conduc¬ 
tivities  of  the  different  rock  layers  are  assumed  to  play  a 
decisive  part,  since  the  wav"  length  in  the  medium  is  more 
reduced  at  higher  conductivity.  The  ratio  between  wave  length 
and  dimensions  will  also  be  decisive. 

(2)  The  measurements  suggest  that  the  antenna  alignment 
with  respect  to  the  coal  seam  and  to  the  transmitter  -  receiver 
direction  is  important.  This  means  that  optimum  propagation 

is  reached  in  a  given  direction  for  a  certain  angle  d’'  which 
must  bo  other  than  0°  or  90*"  • 

Some  mea3\;rement3  have  shown  that  damping  in  the  adjacent  rock 
is  e:.tremely  high  so  that  signals  above  ground  are  practically 
screened  off  atdepths  of  no  more  than  25O  m. 

The  measured  values  of  field  strength  for  distances  of  225, 

X'7C  nir*A  m  ov*A  r--?  A  r\  YTV  nnH  Pin*.  1.^0.  Tonoai'iT'i/r 

the  measurement  at  =  0°  an  i  -y -■  Q3o,  the  value  0,22  was 
determined  f'^r  the  function  p(r,/l}  lefined  in  18].  This 
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corresponds  to  an  average  conductivity  of  o  =  10  Ohm  m 

at  a  frequency  of  3  kc/sec.  Considering  the  highly  heterogeneous 

medium,  this  value  is  to  be  looked  upon  as  a  pure  mean  value. 

This  can  also  be  seen  from  the  values  of  d-c  measurements 
determined  at  several  limited  points. 

1,9  St.  Gertraudi  mine 

Besides  our  measurements  in  German  mines,  we  also  conducted 
detailed  studies  in  nearby  Tyrolese  mines.  The  geological 
conditions  are  well  known  (as  far  as  they  influence  the 
propagation  of  7LF  waves),  and  we  therefore  concentrated  on 
measuring  over  very  large  distances. 

For  this  purpose,  the  1  kw  transmitter  mentioned  above  was  set 
up.  Being  a  mains-supplied  instrument  with  a  power  consumption 
of  almost  3  kw  when  in  full  operation,  it  had  to  be  set  up 
permanently  in  the  mine.  The  position  of  the  laboratory  was 
chosen  so  that  the  distances  accessible  for  measurements  in  the 
St.  Gertraudi  mine  ai'e  as  large  as  possible.  Pig.  1,21  show's 
that  the  maximum  distance  is  approximately  1.5  km.  The  region 
is  especially  suitable,  because  it  contains  no  rails  or  other 
equipment . 

The  advantage  of  a  powerful  but  stationary  transmitter,  however, 
also  entailed  the  drawback  of  the  transmitter  antenna  having  to 
be  fixed  in  a  system  of  galleries  in  order  to  utilize  the  whole 
output  power.  The  antenna  thus  cannot  be  rotated.  The  fixed 
SAVIII  transmitter  antenna  has  been  used  for  most  of  the  measure¬ 
ments  in  St.  Gertraudi.  The  problems  discussed  in  part  III  made 
measurements  over  distances  larger  than  I5OO  m  necessary.  For 
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this  purpose,  a  number  of  measurements  had  to  be  made  at  ground 
level,  i»e,  at  zero  depth.  How  far  these  measurements  are 
meaningful  will  also  be  discussed  in  chapter  III. 

Let  us  now  discuss  a  measurement  made  at  a  frequency  v'  =  Jkc/sec 
(as  In  the  last  annual  report).  The  range  of  distances  was 
approximately  100  m  -  3.3  km  (Fig.  1,2l),  the  measurements  on  the 
western  drift  being  continued  on  the  earth's  surface  from  a 
distance  of  I5  km  onwards.  The  SAVIII  fixed  frame  was  used  as 
transmitter  antenna.  The  accuracy  of  measurement  was  considerably 
increased  by  repeating  the  mea'^urement  and  by  using  several 
receiving  antennas.  The  antennas  EA  IV,  EA  VI,  and  the  calibra¬ 
tion  coil  ESP  described  in  [1j  were  used  as  receiving  antennas. 

The  received  voltage  had  to  be  reduced  to  the  calibration  coil 
in  order  to  make  comparisons  possible,  because  the  three  antennas 
have  different  sensitivities.  These  values  are  given  in 
Table  AV  column  U  and  clotted  in  Fig.  3»1*  \Vhen  measuring  with 
the  fixed  frame,  a  correction  of  the  *!)•  values,  i.e.  reductioi 
of  the  measured  values  to  the  angle  “S'  =  0°  proved  impossible  and 
the  measured  points  plotted  in  Fig.  5»1  ahow  certain  deviations. 
The  present  transmitter  antenna,  which  cannot  make  full  use  of 
the  new  amplifier,  cannot  be  used  for  measuremerta  over 
distances  greater  than  A  km.  There  are  therefore  plans  to  build 
a  much  larger  frame  in  the  near  future  which  will  make  possible 
measurements  over  a  distance  at  least  three  times  as  large 
It  will  then  be  possible  to  conduct  transmission  and  reception 
experiments  between  the  Schwaz  and  St.  Gertraudi  mines,  which 
are  10  km  apart.  Since  the  sensitivity  of  the  receiving  antenna 
may  also  be  improved  (very  promising  experiments  have  already 
been  made),  v/o  may  expect  to  reach  much  greater  iistanceo  still. 


FIG.  1 21 
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Th?  favorable  position  of  the  mines  in  the  Tyrol,  the  maximum 
distances  between  the  individual  mines  being  10  km,  proved  very 
advantageous . 


1-22 


1,10  Determination  of  the  electrical  conductivity  of  rock 

A  four-electrode  configuration  (usually  the  Wenner  configuration) 
was  used  for  determining  the  conductivity  of  various  parts  of 
rock  (mainly  in  mines).  Only  direct  current  was  used  for  this 
purpose.  It  was  the  purpose  of  this  examination  to  determine 
the  conductivity  of  those  regions  in  which  the  propagation  of 
VLP  raves  was  measured. 

1.  Theoretical  considerations: 


According  to  the  analogy  between  an  electrostatic  field  and 
a  steady  electromagnetic  field,  a  point  electrode  in  a 
homogeneous  body  with  the  conductivity  o  produces  a  potential 
given  by  (p  «  in  the  receiving  point  P,  If  there  is  a 

second  electrode  with  the  current  (  -  I )  in  this  medium,  the 
potential  difference  between  the  points  and  (Fig.  1,22) 
is  calculated  as  follows: 
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If  the  electrodes  are  not  in  a  full  space,  but  on  the  surface 
of  an  unbounded  homogeneous  somiopace,  equation  (l)  goes  over 


into 
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or 


U  =  •  k 

2  no 


0  =  -  .  k 
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(2) 


If  the  electrodes  and  ^.s  well  as  the  probes  and 
are  arranged  alv  ’.g  a  straight  line  at  the  distance  a  (V'enner 
configuration),  k  becomes  l/a. 


2.  Measuring  arrangement  (Pig.  1,23  and  Fig.  V) 


External  circuit: 

Two  6-v  lead  accumulators  connected  in  series  feed  a  motor- 
generator  whose  secondary  circuit  provides  a  smooth  direct 
voltage  of  250  V,  This  voltage  is  applied  to  the  electrodes 
and  via  an  ammeter  having  several  ranges  of  measurement.  If 
the  conductivity  of  the  rock  to  be  studied  is  high,  the  12  v  of 
the  accumulators  instead  of  the  25O  v  generator  voltage  can  be 
applied  to  the  electrodes  via  the  instrument  by  means  of  a 
change-over  switch. 


Voltage  measurement 

The  voltage  U  between  the  probes  3^  and  has  to  be  measurea 
with  a  compensation  circuit  in  order  to  eliminate  the  high 
contact  resistance  (50  -  300  kohms)  of  the  electrodes  and  probes. 
A  stabilizing  circuit  equipped  with  two  transistors  AC  125  nnd 
a  Zener  diode  OAZ  203  provides  a  constant  voltage  of  10  or  1  v 
across  a  50  kohms  helical  potentiometer  (type  Spektrol  860). 

The  compensation  voltage  on  the  slider  of  the  potentiometer  is 
lead  to  the  probe  via  a  100  pA  zero  indicator  connected  with  a 
transistorized  d-c  amplifier,  before  the  input  of  th«  amplifier 
there  is  a  second  compens?- 1 i^n  circuit  for  comj.ensating  the 
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contact  voltage  and  polarization  voltage  of  the  probes  with  the 
electrode  current  being  switched  off.  The  2  kohms  potentiometer 
with  central  tapping  has  a  logarithmic  characteristic  so  that 
the  small  values  occurring  frequently  can  be  adjusted  accurately. 

The  amplification  factor  of  the  d-c  amplifier  is  20.  Thus,  a 
■■8 

value  of  5*‘*0”  A  can  easily  be  read  on  the  instrument.  If  an 

extreme  resistance  of  1  Mohm  is  assumed  between  the  probes,  a 

_2 

voltage  change  of  5*10  v  is  necessary  to  produce  a  sufficient 
deflection  of  the  zero  instrument.,  Frequently,  however,  an 
accuracy  in  voltage  measurement  of  5*10  ^  v  can  be  reached  owing 
to  thf  use  of  a  built-in  key  and  owing  to  the  frequently  small 
resis  :ance  between  the  probes. 


3.  Measuring  accuracy: 

If  the  quantities  J,  U,  r^ ,  and  rg  have  errors  p^,  Pj^ 
the  error  of  o  will  be 


If  the  errors  are  substituted  in  (5)  and  if  the  measuring  error 
of  the  unit  distance  p^  is  substituted  for  the  error  p^  of  the 
distance  measured  according  to  p^-p^'^r,  we  obtain 


1  cf.  2  c^I^  2  2  /J_  T  T  _1_  \ 

'12^1'*’  4^'*'  22  ^e  '3*^  3'*’  3'*’  3  ) 

i(u^  4n''u  ®  r^  v'  R? 


The  error  for  the  Wenner  configuration  will  thus  be 


^  =  ^\\T!  Mr;  4'I 


(4) 


(R^  *  r2  »  a,  r^  -  R^  =  2a). 
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For  the  conducted  measurements,  the  following  maximum  values  are 
to  be  substituted: 

-f  =  0,02  =  0,05  =  0,02  a  =  0,5 

V-^ith  these  values,  the  relative  error  is  0,07.  Thus,  the 
measurements  in  the  most  unfavoiable  case  will  have  an  error 
of  7^-. 

4.  Measurements: 


Electrodes  and  probes : 

For  electrodes  as  well  as  for  probes,  holes  were  made  into  the 
rock  by  means  of  a  spudding  bit  fed  by  an  accumulator,  which 
filled  with  cotton  wool  soaked  in  a  salt  solution.  Then, 
steel  studs  were  driven  into  the  cotton  wool  to  fasten  the 
wires  for  the  measuring  device.  Although  this  method  of 
contacting  yield's  very  low  contact  resistances,  the  production 
of  electrodes  for  a  large  number  of  measuring  points  would  have 
consumed  too  much  time.  Commercial  bolt  setting  guns  (Hilti 
DX  100)  were  used  to  set  steel  nails  into  the  rock  for  the  below 
measurements.  This,  however,  was  only  possible  in  small  cracks 
in  the  rock.  Two  suitable  sites  were  chosen  for  installing  two 
Wenner  configurations  at  a  vertical  distance  of  approximately 
10  cm  in  order  to  study  whether  the  results  of  the  measurements 
are  affected  by  a  possible  difference  in  conductivity  in  the 
cracks  filled  with  wet  dust  (Fig.  1,24)» 

abed 

+  +  +  + 


a  *+ 
C 

Fig.  1  ,24 
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The  studs  a,  h,  Cf  and  d  were  placed  into  cracks,  whereas  the 
atuds  A,  B,  C,  and  D  were  set  into  solid  rock,  their  combination 
being  varied  during  the  measurement.  The  effect  of  the  cracks 
carrying  the  electrodes  was  assumed  to  be  negligible,  as  the 
deviation  from  the  mean  value  is  < 

The  first  measurement  was  made  to  determine  the  conductivity  over 
a  large  distance  by  means  of  the  four-electrode  method.  For  this 
purpose,  the  configuration  was  chosen  such  that  it  seemed 
justified  that  the  expression  for  the  full  space  should  be 
applied.  A  map  is  given  in  Fig.  1,25. 

Measurements  were  conducted  with  probes  on  the  ceiling  and  on 
the  walla  of  the  gallery.  The  electrodes  were  made  different, 
measurements  on  these  electrodes  were  repeated  several  times 
after  periods  of  some  weeks,  as  the  water  level  in  the  gaps  of 
this  mine  changes  considerably.  The  below  mean  value  of  52 
measurements  was  obtained  for  the  apparent  specific  conductivity 
on  this  measuring  line. 

Og  *  1.62»10~^  mho/m 

the  deviation  from  this  mean  value  being  CBffo,  It  has  still 
to  be  found  out  (by  changing  the  distances  of  the  configuration) 
whether  this  value  represents  the  actual  d-c  conductivity.  For 
this  purpose,  similar  measurements  over  smaller  distances  were 
also  made. 

First,  studies  were  necessary  to  find  out  whether  the  measuring 
results  are  affected  by'  a  difference  in  conductivity  of  the 
often  thick  and  damp  layer  of  dust  and  the  rock  below.  For  this 
purpose,  the  theory  of  the  stratified  medium  with  a  simplifica¬ 
tion  for  only  one  layer  was  used.  The  results  are  shown  in  Fig. 


1 ,26. 
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The  symbols  have  the  following  meaning: 


0-0, 


0  +  0, 


a 

^  h 


p  =  resistivity  of  covering  layer 

»  resistivity  of  bottom 

p  »  apparent  resistivity  measured  * 

■;  cm  was  used  as  the  smallest  value  of  a.  Studs  mounted  to  a 
pertinax  rod  and  pressed  toward  the  rock  by  springs  were  used 
as  electrodes  and  probes.  Fig.  1,2?  gives  the  result  of  one 
measurement.  The  ordinate,  however,  gives  p  instead  of  v  and 
the  abscissa  gives  a  instead  of  u. 

Hence,  a  thin  surface  layer  actually  has  a  higher  conductivity, 
which  is  ineffective  already  at  a  =  0.30  m.  \7e  only  conducted 
measurements  with  larger  a. 

The  neighborhood  of  the  electrodes  S^,  and  £2  studied  in 
detail.  Fig.  1,28  shows  a  diagram  of  the  used  electrodes. 

Various  combinations  of  the  electrodes  were  used  as  four-electrode 
configurations,  and  the  corresponding  apparent  conductivity  was 
determined.  Table  XVI  shows  the  results; 


These  values  differ  considerably  from  each  other,  the  smallest 
value  being  0.74 mho/ra  and  the  largest  value  leing 
5.00*10”^  raho/m.  The  derived  error  of  measurement  was  checked 
by  repeated  measurements  showing  the  results  to  be  well 
reprod-v-ible.  For  equal  electrodes,  however,  the  values  on  the 
whole  increase  if  the  probes  are  placed  subsequently  at  FG,  PH, 
and  HI.  Owing  to  this  regularity,  the  apparent  conductivity  is 
assumed  to  have  different  values  for  two  reasons: 

(1)  There  are  differences  in  the  actual  specific  conduc¬ 


tivity  of  the  rock. 
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(2)  The  differences  r.re  due  to  the  effect  of  the  gallery  . 
which  with  a  cross  section  of  approximately  2  m  by  2m  is  equally 
large  as  the  distances  of  the  used  configuration. 

Some  other  measurements  prove  the  two  cor.iclusions  to  be  correct. 
Four  holes  at  a  distance  of  two  meters  each  (depth  2  m)  were 
drilled  into  the  wall  of  a  larger  mine.  The  V/enner  config^xration 
was  used  for  determining  o^  at  different  distances  of  the 
electrodes  and  probes  from  the  surface.  If  o^  is  always  calculated 
with  the  expression  of  a  semispace  regardless  of  the  depth,  it  has 
to  increase  with  the  depth  provided  the  rock  is  homogeneous.  The 
result  (given  in  Table  XTIl),  however,  at  this  point  is 
independent  of  the  depth,  i.e.,  the  actual  conductivity  increases 
with  the  depth.  This  is  confirmed  by  a  measurement  in  which 
Wenner  configurations  above  this  point  were  shifted,  with  the 
distance  a  being  constant  (a  =  0.5,  1,  1.5»  2,  and  2.5  m).  The 
expected  increase  in  a  with  a  is  shown  in  Pig.  1,29.  Every 
value  is  plotted  in  the  diagram  at  the  point  M  which  is  the 
center  of  the  configuration.  Hence,  the  conductivity  a  changes 
irregularly  within  the  above  limits  already  at  small  distances, 
vertical  to  the  configuration  as  well  as  along  it.  Similar 
results  were  found  also  at  other  measuring  sites. 

The  measurement  below  shows  the,t  the  effect  of  the  gallery  cannot 

always  be  neglected.  Each  of  the  two  side  walla  of  a  gallery 

was  provided  with  a  V'enner  configuration  having  a  =  1.5  m.  The 

site  wes  chosen  such  th/.t  was  the  same  on  either  side  (Fig.  1,50 

The  value  obtained  for  a  is  different,  if  the  probes  and 

s 

electrodes  are  placed  at  different  lides  of  the  gallery. 
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calculated,  from  the  expression  (1). 


— - - - - 1 - 4- 

a  h  c  d 


a*  b'  C  d* 

Fig.  1,30 

The  analogous  values  for  &  =  ^  m  coincide  much  better.  Because 
of  the  inhomogeneity  of  rock,  this  measurement  was  conducted  for 
several  points  by  shifting  the  configuration  always  by  5  w* 

Fig.  1,31  shows  the  result.  The  difference  In  the  deviation  of 
the  corresponding  points  of  measurement  is  probably  due  to  an 
asymmetry  in  homogeneity  with  respect  to  the  gallery.  The 
values  measured  at  the  same  site  with  a  =  10  ra  are  in  sufficient 
agreement  with  allowance  for  the  error  limits  plotted  in  one  of 
the  curves  (Fig.  1,52).  This  is  assumed  to  show  that,  under 
the  given  conditions,  the  effect  of  the  space  of  the  gallery 
from  this  a-value  onward  can  be  neglected  and  that  the 
expression  (l)  yields  the  correct  values  for  the  conductivity 
in  homogeneous  rock.  In  order  to  get  a  better  idea  of  the 
restrictions  that  are  to  be  observed  for  determining  the 
conductivity  in  mines  by  the  '•"enner  configuration  it  is 
desirable,  that  a  general  relation  among  the  three  parameters 
cross  section  of  the  galleries,  a  and  o  be  derived.  This  is 
planned  to  be  done  in  the  continuation  of  this  study. 

After  measuring  in  various  parts  of  the  mine  with  allowance  for 
the  above  experience,  we  may  now  say  that  the  d-c  conductivity 
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changes  within  the  values  0.5*  10"^  and  5 •0*10'’^  mho/m  given  above, 
the  homogeneous  ranges  being  but  small. 

5*  Furthermore,  the  results  of  measurements  in  German  and  Dutch 
mines  which  we  visited  in  summer  1965,  shall  be  summarized. 

( ' )  Salzgitter,  mine  Konrad 

A  highly  homogeneous  ore  body  v/hose  conductivity  is 
5.1 *10”^  mho/m.  At  the  edge  of  this  body  where  it  is  mixed 
with  dead  rock,  the  conductivity  decreases  to  1.2*10”^  mho/m, 
in  dead  rock  (in  the  1200  m  gallery)  to  3*0*10*’  raho/m. 

(2)  Hannover,  potash  mine  Hansa  III-Empelde 

Measurements  with  the  available  device  were  impossible  because 
of  the  high  resistivity.  In  any  case,  the  conductivity  is 
smaller  than  1*10“^  mho/m, 

(5)  Siegerland,  ore  mines  Fuesseberg  and  Willroth 
The  ore  of  both  mines  is  highly  inhomogeneous  and  mixed  v;ith 
dead  rock.  The  results  in  ore  and  also  in  dead  rock  lie  between 
2»10"’^  mho/m  and  8»10"^  mho/m. 


(4)  Heerlen,  coal  mine  Oranje  -  Nassau  III 

Coal  showed  a  clear  diffei’ence  between  the  results  in  the 

direction  of  the  seam  and  those  vertical  to  the  seam. 

Result;  in  direction  of  the  seam  a  =  8.0*10"  mho/m 

vertical  to  the  seam  o  =  1.S»10  ^  mho/m 

The  possibility  of  measuring  in  dead  rock  was  restricted  because 

-3 

of  the  iron  props;  the  few  measured  values  lie  between  3*3*10 

•2 

and  1,1 *10”  mho/m. 
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6,  Specimen  measurements,  frequency  dependence. 

So  far,  the  frequency  dependence  has  only  been  measured  for 
specimens.  For  making  the  specimens,  a  cylinder  with  a 
diameter  of  5  cni  was  bored  out  of  a  larger  piece  of  rock  (Fig.  VI ) 

and  was  then  cut  into  disks  0.2  -  1 .0  cm  thick  with  a  stone 

saw  (Fig.  VIl),  The  loss  factor  and  the  capacitance  of  the 
samples  were  measured  in  a  measuring  capacitor  (Fig.  VIII )  built 
for  this  purpose  by  a  capacitance  bridge  (General  Radio, 

■type  716  -  C)  under  different  conditions.  From  these  values, 
a  and  £  were  calculated.  The  humidity  and  the  porosity  of  the 
specimens  were  determined  from  the  difference  in  weight  before 
and  after  drying  at  125°C  for  24  hours. 

In  the  present  report  we  only  want  to  discuss  a  few  results  in 

brief,  namely  the  values  obtained  for  two  specimens  (D^  and  Pg) 

from  the  St.  Gertraudi  mine,  two  ore  specimens  (S^  and  S2)  from 
the  Salzgitter  ore  mine,  and  tv/o  specimens  (G^  and  G2)  from  the 
V/illroth  mine.  The  two  specimens  are  always  chosen  from  a  large 
number  of  specimens  being  nearly  equal  with  respect  to  their 
electrical  properties.  and  D2  consist  of  dolomite  of  low 
porosity  (pore  volume  approximately  O.IZfo  of  the  total  volume), 
and  S2  are  of  oolitvic  ore  of  much  greater  porosity  (pore 
volume  approximately  14f ),  and  G^  and  G2  are  ore  specimens  mixed 
with  dead  rock,  whose  porosity  is  somewhat  higher  than  that  of 
dolomite  differing  considerably  in  various  specimens  -  in  contrast 
to  the  specimens  D  and  S.  Thve  electrical  properties  7/=re  first 
measured  immediately  after  cutting  the  specimens,  i.e.  in  a 
largely  natural  state.  The  specimens  were  then  dried  and  the 
measurements  were  repeated.  The  third  measurements  were  made 


with  the  specimens  saturated  in  distilled  water.  The  val  iS  of 
t  and  a  are  given  in  Table  XIX  and  Pigs.  1,35  -f-nd, ■’!,34 

The  values  for  a  are  all  several  powers  of  ten  smaller  tha^n 
those  measured  with  direct  current  in  outcrops.  Thus  it  seems 
that  boring  and  cutting  the  specimens  changes  their  electrical 
properties  considerably,  and  henco  that  t!^e  frequency  curves  for 
the  conductivity  measured  for  the  specimens  are  completely 
different  from  the  d-c  values  measured  in  outcrops.  On  tho 
other  hand,  since  the  factor  between  the  values  at  100  cps  and 
10  kc/sec  is  two  powers  of  ten  at  least,  it  seems  necessary  that 
the  frequency  dependence  should  also  be  measured  at  the  outcrop. 
For  this  purpose,  a  second  transmitter  was  set  up  in  the 
St,  Gertraudi  mine  lab  (besides  the  Savage  amplifier  used  for 
propagation  experiments)  v^hich  also  had  an  output  power  of  1  kw 
and  was  suitable  as  a  voltage  source.  However,  our  work  on  the 
measurement  of  the  frequency  dependence  of  conductivity  in  the 
outcrop  has  not  progressed  far  enough  for  the  results  to  be 
discussed  in  the  present  report. 


TABLE  I 


Calcul2,ticn  of  the  function 


f  =  5  kc/sec 


TABLE  II 


Function 


f  (r,A) 


H 

10“^ 

10“^ 

2‘ }0"^ 

l5*io"^ 

1 

10 

2.02 

1.95 

1.83 

1.47 

1.16 

20 

1.97 

1.44 

1.26 

0.88 

0.72 

50 

1.71 

0.75 

0.54 

0.5? 

0.261 

100 

1.16 

0.37 

0.29 

0,182 

0.12: 

200 

0.65 

0.2 

0.145 

0.09 

- 

500 

■  ■  - - i 

0.26 

0.081 

0.059 

_ _ 1 

- 

- 

TABLE  III 


Measiiremant  I 

^^=0°  50-' 


t 

U 

e 

Ue 

‘  y 

B 

fc 

© 

0° 

5.3 

5.3 

5.57 

5 

5 

5.16 

K 

4.5 

4.6 

4.65 

4.8 

4.5 

4.46 

60° 

2.5 

2.8 

2.68 

9^0 

2.5 

3.2 

2.58 

90° 

0.85 

0.9 

0 

120 

1.7 

1.7 

0 

120 

3.0 

2.6 

2.68 

150° 

2.5 

5.2 

2.58 

15G° 

4.8 

4.6 

4.65 

180° 

^t.3 

4.5 

4.46 

180 

5.4 

5.3 

5.37 

210° 

5.2 

5.16 

350° 

4.8 

4.6 

4.65 

0° 

4.2 

4.3 

4.46 

300° 

5.0 

2.8 

2.68 

350° 

2.4 

3.2 

2.58 

270 

1.0 

0.9 

0 

2.5 

2.8 

2.68 

210° 

4.5 

4.6 

4.65 

0° 

5.3 

5.5 

5.37  } 

u 

Q 

(;o“ 

“e 

K 

r 

'^=  90° 

U 

e 

u 

e 

! 

1 

60° 

4.7 

4.7 

4.72 

270° 

4.6 

4.6 

4.75 

4.2 

4.1 

4.16 

60° 

4.3 

4.3 

4.12 

0° 

2.7 

2.5 

2.41 

5^0 

2.7 

2.5 

2. 38 

330° 

0.8 

0.8 

0 

4.6 

4 .6 

4.75 

300° 

2.1 

2.5 

2.41 

120 

4.0 

4.3 

4.12 

270° 

4.0 

4.1 

4.16 

150° 

2.0 

2.5 

2.38 

qo° 

4.0 

4.1 

4.16 

180 

0.6''^ 

0.6 

n 

120° 

2.4 

2.5 

2.41 

210° 

2.7 

2.5 

2.38 

150° 

0.8 

0.8 

0 

240 

^..3 

4.5 

4.12 

180 

2.6 

2.5 

2.41 

TABLE  IV 


Measurement  II 
r  =  25  m 


y 

U 

e 

0 

0 

U 

e 

O' 

i 

30° 

5. 

5e 

0° 

60° 

90° 

120° 

150° 

180° 

330° 

300 

270° 

2.4 

2.0 

1.3 
0.85 
1.6 

2.2 

2.4 

2.2 

1.55 

0.85 

2.4 

2.2 

1.6 

0.85 

1.6 

2,2 

2.4 

2.2 

1.6 

0.85 

2.5 

2.16 

1.25 

0 

1.25 

2.16 

2.5 

2.16 

1.25 

0 

0° 

60° 

90° 

120 

150° 

180° 

300° 

2.2 

2.4 

2.3 

1.8 

1.5 

1.7 

2.2 

1.7 

1.5 

2.2 

2.4 
2.2 
1-7 

1.5 
1.7 
2.2 
1*7 
1-5 

2.34 

2.7 

2.34 

1.35 

0 

1.55 

2.34 

1.35 

0 

6u 

90° 

r 

U 

e 

tj 

e 

®e 

y 

U 

Q 

®e 

lo 

60° 

5.4 

5.4 

3.5 

90° 

5.2 

3*2 

5. 18 

90 

2,8 

2.8 

2.85 

60° 

2.7 

2.7 

2.74 

120 
^  <-^o 

1.7 

1.7 

1 .65 

1.75 

1.6 

1.59 

150 

0.55 

0.55 

0 

0° 

1.55 

1*5 

0 

180 

nO 

1 .6 

1.7 

1.65 

120° 

2.7 

2.7 

2.74 

210 

■rnO 

2.8 

2.8 

2.85 

150° 

1.6 

1.6 

1.59 

30 

2.8 

2.8 

2.85 

180° 

1.3 

1*5 

0 

0 

1.75 

1.7 

1.05 

210° 

n 

1.6 

1.6 

1.59 

2.,0 

2.7 

2.7 

2.74 

Table  V 

transmitter  at  PPIO7 
f  »  3  kHz;  I  =  0.75  a;  SAVI 


pp 

m 

'e  W 

51 

455 

15° 

1.5 

64 

558 

42° 

1.5 

99 

630 

400 

1.1 

10 

420 

52° 

1.5 

52 

290 

150 

1.5 

95 

240 

0° 

43 

100 

627 

6° 

1.0 

105 

300 

10° 

1*5 

Table  VII 


measurement  1 ;  transmitter  PP9 
f  =  5  kc;  X  =  3.5  a;  SAVI; 

receiving  antenna  EAVI 


PP 

m 

a 

8 

67 

40° 

5500 

7 

140 

250 

200 

6 

215 

150 

55 

5 

285 

5° 

27 

5.1 

262 

00 

26.5 

3.2 

460 

10<^ 

4.5 

2 

535 

10° 

2.8 

2.1 

520 

0° 

2.4 

Table  VI 

transmitter  at  PPIO7S  f 


4° 

^eLl'Xl 

0 

4 

4.10 

50 

3.6 

3.54 

60 

1.9 

2.05 

90 

0.3 

0 

120 

2.1 

2.05 

■150 

3.6 

5.54 

350 

3.6 

5.54 

300 

2.0 

2.05 

2V0 

0.3 

0 

kc;  r 

=  240  m , 

receiver  PP95 

60° 

1 0 

4 

I'eH 

n.UY! 

270 

2.0 

1.92 

500 

1.4 

1.11 

550 

0.4 

0 

0 

1.0 

1.11 

50 

1.9 

1.92 

60 

2.3 

2.22 

240 

2.2 

2.22 

210 

1.7 

1.92 

340 

0.4 

0 

Table  VIII 

measurement  1 ;  transmitter 

f  =  3  kc;  I  =  0<4  a,  SAVI 

PP  m 

Table  X 

measurement  2:  transmitter  Sg 
f  =  3  kc;  I  =  0.7  a,  SAVI 

PP  m  ip,fuVj 

1 

55 

52.7 

505 

60  57 

210 

605.1 

80 

16.5 

c 

160  58 

60 

4.1 

150 

10 

603.1 

164  46 

53 

4 

155 

9.2 

d 

222  19 

26 

^1 

156  22 

42 

1 

125  28 

60 

4.1 

81  57 

150 

Table  IX 

measurement  2s 

transmitter  S  j  f  =  5  kc 

;  SAVI;  I  = 

5.5  a; 

receiving 

antennas  BAIV 

and  EAVI 

PP 

m 

U^(EAVI) 

Ug(EAIV) 

1 

170 

45 

6 

50 

6,05 

1.1 

119 

42 

11 

45 

10.1 

2 

141 

55 

11 

55 

10.9 

2.1 

166 

52 

15 

67 

14.5 

3.2 

122 

68 

70 

450 

81 

4 

107 

86 

100 

550 

106 

5 

110 

71 

57 

170 

56 

5.2 

109 

67 

20 

100 

20.2 

6 

126 

54 

17 

80 

16.7 

7.1 

15? 

42 

11 

52 

10.8 

7 

150 

45 

15 

55 

12.^ 

8 

180 

54 

8 

- 

8 

8.1 

18  2 

54 

7 

55 

7.1 

9 

202 

50 

4 .6 

25 

4.9 

Table  XI 


o« 

r  = 

262  m 

90° 

4,0 

4*° 

0 

24 

24.1 

265 

18 

17.8 

350 

21.5 

20.8 

240 

16.5 

16.1 

300 

14 

12 

210 

11.5 

10.4 

270 

0.5 

0 

175 

0 

0 

240 

10.5 

12 

180 

5.3 

1.55 

210 

10 

12 

500 

14.2 

14.6 

180 

24 

24.1 

330 

7 

7.52 

30 

20 

20.8 

0 

3 

1.55 

60 

11 

12 

3Q 

11 

10.4 

90 

0.5 

0 

60 

16.5 

16.1 

120 

14 

12 

90 

17.5 

17.6 

150 

21.5 

20.8 

150 

6»4 

7.52 

Table  XII 

transmitter 

PPS^,  r 

“  55  ni,  f 

=  5kc 

5  -  0° 

90° 

4,0 

4,0 

"eM 

100 

53 

52.7 

270 

50 

30 

150 

47 

45.5 

240 

27 

26 

120 

20 

26.4 

210 

16 

15 

90 

1 

0 

180 

2 

0 

60 

23.5 

26.4 

150 

14 

15 

30 

45 

45.5 

120 

25 

26 

210 

45 

45.5 

90 

30 

30 

240 

24 

26.4 

270 

1 

0 

0 

0 

27 

26.4 

Table  XIII 

r  a  150  m,  transmitter  PPS^ 


"eH 

4" 

U  '^V^ 

e  f*  - 

210 

10.5 

10.15 

270 

10 

10 

180 

9 

8.8 

240 

8 , 6 

8.65 

150 

4.5 

5.07 

210 

5.0 

5.0 

120 

0.2 

0 

180 

0.8 

0 

90 

4.8 

5.07 

150 

4.6 

5.0 

240 

8.8 

8.8 

120 

8.4 

8.65 

270 

5.0 

5.07 

90 

9.6 

10 

300 

0.2 

0 

300 

8 . 6 

3. 65 

Table 

XIV 

Table  XV 

transmitter: 

f  = 

3  kc, 

transiiiitter :  f  = 

3kc,  I  =  25a,  SAVI 

I  a 

la-,  SAVI 

m 

+  10/. 

r 

r 

185  39 

72.5 

280  33 

21  .84 

225 

0 

1 

380  21 

9.3 

225 

90 

0.5 

475  50 

4  =  94 

375 

0 

0.5 

580  30 

2.7.4 

375 

90 

2.3 

630  50 

1 .71 

500 

0 

0.2 

720  50 

1 .215 

040  50 

0.85 

840  50 

0.71 

960  50 

0.545 

1040  30 

0.45 

1100  23 

0.445 

1200  28 

0.35 

1270  27 

0,344 

1320  27 

0.262 

1550  23 

0.543 

1C80 

0 

0,151 

2000  27 

0.0813 

2200  24 

0.1?7 

3000  3 

■6 

0.058 

3300 

4 

0.0197 

Table  XVI 


electrodes 

probes 

k 

J 

IT 

cj’IO”^  mho/i 

B  -  K 

P  -  G 

0.246 

1.75 

0.17 

2.02 

P  -  H 

C.4C5 

1 .60 

0.25 

2.17 

P  -  I 

0.551 

1.55 

0.21 

2.06 

G  -  I 

0.1C5 

1.55 

0.06 

2.16 

G  -  H 

0.179 

1 .60 

0.10 

2.29 

I  -  H 

0.074 

1.45 

0.05 

1 .71 

F  -  K 

G  -  H 

0.275 

1.15 

0.12 

2.10 

G  -  I 

0.222 

1 .20 

0.11 

1.93 

I  -  H 

0.053 

1.25 

0.032 

1 .65 

a  -  K 

H  -  I 

0.06^ 

1.10 

0.075 

0.74 

D  -  H 

P  -  G 

0.379 

5.40 

0.60 

1 .71 

E  -  I 

P  -  G 

0.255 

6 .60 

0.67 

1 .99 

P  -  H 

0.556 

6.70 

0.86 

2.21 

G  -  H 

0.103 

6.60 

0.18 

5.00 

E  -  L 

P  -  G 

0.232 

3.50 

0.32 

2.02 

F  -  H 

0.550 

5.70 

0.44 

2.55 

P  -  I 

0,425 

3.70 

0.52 

2.40 

H  -  I 

0.073 

5.60 

0.06 

2.62 

G  -  I 

C.ini 

5.65 

0.19 

2.92 

j’  -  :• 

0  -  ^ 

0.214 

2.60 

0.20 

2.21 

G  -  I 

O.3C8 

2,60 

0.28 

2.28 

H  -  I 

0.094 

2.60 

0 .09 

2.16 

B  -  I 

0.  .62 

2.40 

0.14 

2.21 

G  -  L 

H  -  I 

0.210 

2.60 

0.29 

T.50 

B  -  1 

P  I 

0.176 

3-50 

0.225 

2.18 

A  -  I 

0.373 

5«40 

0.43 

2.35 

f 


Table  XVII 


Depth  of 

0 

0.5  1.0 

1.5 

2.0 

configuration 

a  *10”^ 
s 

2 

.02 

2.12  2.27 

2.05 

2.01 

Table  XVIII 

Electrodes 

Probes 

k 

J  U 

10"^ 

a  -  d 

b  -  c 

0.667 

1.62  0.321 

2 

,68 

a'-d' 

b'-c' 

0.667 

6.15  1.245 

2.62 

a  -  d 

b '  "-c ' 

0.349 

1.62  0.050 

9 

.03 

a'-d' 

b  -  c 

0.549 

6.15  0.192 

8 

.87 

Table  XIX 

A.  Dried  specimens 

1) 

D 

^  (Fig* 

1,33;  D^I) 

Dg  (Pig. 

1,33;  3)2^) 

f [kcpsj 

L 

r 

t 

S' 

10 

4.11 

9.25*10‘^ 

4.92 

1.25*10"® 

100 

4.09 

1.22.10”^ 

4.88 

1 .42.10"'^ 

2) 

S 

1  (Fig. 

1,34;  S^l) 

Sg  (Fig. 

1,34;  s^l) 

f [kcpsj 

6 

r 

t 

6" 

0.1 

23.7 

8.4.10"^ 

50.0 

9.70.10"^ 

0.2 

22.1 

1.83.10"® 

49.7 

2.16.10"® 

0.7 

20.2 

6.90.10"® 

48.8 

6.06.10"® 

1 

19.7 

1 .00.10"^ 

48.6 

1 .01.10"”^ 

2 

18. 7 

1.87.10"^ 

48.0 

2.60*10"^ 

5 

17.4 

4.55*10"^ 

46.5 

9.80.10"^ 

n 

1 

^  /X 

1  1  ,v 

5.56.10  ^ 

45.6 

.  . _  . - -6 
1 . 5U. 1u 

10 

16.7 

7.57-10"'^ 

44.7 

2.54-10"® 

IOC  • 

14.8 

5.96.10"® 

58.1 

3.33.io"5 

Table  XIX 


B. 

1) 


Specimens  saturated  with  distilled  water 


f  [kcps] 


Dt,  (Fig.  .1-,55; 


Do  (Fig.  1,53;  D.II) 


0.1 

o 

• 

C 

2.45-10"® 
^  “8 

0.2 

28.4 

5.13*10 

0.5 

25.2 

1.51‘10~’^ 

&-.? 

24.0 

2.15*10"^ 

1 

22.6 

3.06-10"’^ 

2  ■ 

20.1  • 

5.74*1.0"'^ 

5  ■ 

16. 9  ■ 

1.45.1.0"® 

7  ■ 

15.8  ■ 

2.02-10"® 

-6 

10  ’ 

14.4  ■ 

2.88.10 

100 

7.2 

1  .71-10"^ 

65.0 

■  1.23-10"’^ 

53.8 

2.10*10"’^ 

43.8 

4.58.10"'^ 

41.2 

5.81*10"^ 

-38>5 

7..,99-ie’‘'^ 

—6 

33.0 

1.47-10 

27.9 

3.46*10  ® 

22.8 

4. 61 *10"® 

19.4 

6.28*10"® 

7.8 

1  n 

2.42*10"^ 

2)-- 

• 

(Fig. 

.1,34;  s^ii) 

Sg  (Fig. 

1,34;  8211) 

f  [kcps] 

£ 

r* 

6- 

0 

■ 

•  1 

.  92.5 

<4.72*10"® 

<-^.62*10"® 

2 

60,6 

5*53*10"® 

186.0 

1.55*10"^ 

5 

57.9 

6.37*10"®  . 

115.7 

1 .87 *10" 5 

,-7^ 

34.0 

6.85*10"® 

105.5 

2.04.10"^ 

10 

51.0 

7;75*io"® 

94.8 

2.50.10"^ 

.  . 

100_ 

23.4 

1.85.10"® 

69.8 

5.64*10"^ 

■5)  • 

• 

G,  (Fig 

*.  1,33) 

Gg  (Fig. 

1,53) 

• 

f  [ikcps^  , 

"  '-®1  3 

>l«i3  ■ 

* 

0.1  • 

66 . 0 

■  1.12.10“^ 

74.9 

7.65* 10"^ 

0.2  • 

65.4 

2.26*10"® 

70.1 

1.52*10"^ 

0.5  • 

64.1 

8.47*10"® 

62.3 

4.12*10"'^ 

* 

1 

63.2  ' 

2.42*10"^ 

56.7 

8.93.10"^ 

2 

61.8  . 

8.06.10"^ 

/ft  c. 

t.  »-  w" 

f. 

5  ■ 

57.2 

2.81.10"^ 

35.3 

5.02*  10"® 

10 

■  51*9 

7.04.10"® 

26.4 

9_.O8=.10"® 

20 

45.7 

1 .55*10"^ 

19.6 

1 .40.i0"5 

50 

35.7 

3.54.10"^ 

15.4 

2.11 • lo“5 

ICO 

29.2 

7.04.10"^ 

11 ,2 

5.05. io"5 
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II.  Theoretical  part 


2,1  Effect  of  Quadrupole  radiation 


The  following  expression  was  obtained  in  [1  ]  for  the  field 

strength  of  a  quadrupole  (p.  17,  expressions  (56)  -  (59))* 

2 

kl  3ik^ 

.  3a^(3cos‘^T9'-l)  expC-ik^r)  {  -  -5-  +  — ^ 


2  ■  5  ■  4 

r  r 


H.^(r,^)  =  2  ®i^2''^exp(-ikgr)  I - 


^  ^  ± 
,2  "  r  \4 


In  order  to  examine  what  effect  these  terms  have  on  the  dipole 
field,  we  take  the  square  of  the  whole  field 

(h®  +  »  Ih^I^  + 


where  ^  is  the  corrective  function  to  be  estimated,  which 
is  calculated  as  follows: 


2 

The  functions 

for  f  »  3  Wc/sec  and  dt  *= 

2.43*10"^  (see  Table 


were  numerically  computed 
'kfTnTTF'  .-7 

\  =  7.69*10  2.43*10  and 
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Table  2,1 


L\ 

7. 69-10“"^ 

2.43*10"^ 

2. 43*  10”'^ 

r[m] 

2? 

0.20 

0.20 

0.20 

10 

1r 

0.002 

0.002 

0.0059 

100 

2.2*10~^ 

9.2-10"^ 

2.6.10"^ 

1000 

5.24 

3.24 

3.24 

10 

0.05 

0.05 

0.04 

100 

2. 56, -10“^ 

2.25*10"“^ 

7.5*10’^ 

1000 

0.9 

0.9 

0.9 

10 

0.09 

0.09 

0.138 

TOO 

0.009 

0.015 

0.081 

10C0 

-3.6 

-3.6 

-5.6 

10 

-0.36 

-0.45 

-1 .98 

100 

-0.12 

-0.19 

-1.73 

1000 

Substitution  of  the  values  for  r  =  100  m,  a  =  5*10  ^[mhos/m'J 
and  f  »  3  kc/sec  yields 

Q  ^  =  0.57  I  =  O'’ 

Q  ^  =  0.008 |h®|^  'S'  =  90'> 

under  the  arbitrary  assumption  that  a^/m  =  t. 

Owing  to  this  result,  the  quadrupole  terms  are  assumed  to  be 
negligible  freia  a  distance  of  100  m  onward,  since  the  assumption 
a^/m  »  1  certainly  describes  the  maximum  effect  of  the 


quadrupole 
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2,2  Effect  of  cavity 

In  our  last  annual  report  M  j,  the  effect  of  the  cavity  in 
the  form  of  an  influencing  secondary  field  in  the  surrounding 
medium  has  been  dealt  with.  This  theory  which  in  the  considered 
parameter  range  yields  no  measurable  effect,  was  examined 
critically  and  the  result  has  again  been  confirmed.  V/e  now 
want  to  give  a  br  * survey  of  these  detailed  considerations: 


From  the  sol  >n  of  the  wave  equation 


k  +  AA^j  =  0 


we  obtained 


At  =  (a 
Ir  '  o  r 


exp(-ik^r)  sinker 

^  3  -  )  C08V 


k^r 


exp(-ik^^r) 

A,^.  =  C  - cos  T> 

Hr  r 


" 


exp(-ik-r)  sinker 

-(a  - - +  B  -  )  sim^ 

'or  .  ' 

kjr 

exp(-iky^r) 

A_^  =  C  - — —  sini'^ 

II'  r 


where  =  wave  number  of  She  medium  in  the  cavity  and 

kj^  =  wave  number  of  the  surrounding  medium.  The  conditions 

of  continuity  are  used  to  eliminate  B.  For  C  and  A^  an  equation 

is  obtained  which  gives  the  effect  of  the  cavity  in  the  form 

of  the  change  in  C  with  respect  to  A.  Prom  [1 '[  (p,  i9)» 

obtain 

E  =  3in(kjr^)/kjr^  -  cosCk^r^)  k^r^  =  x<^1 


yielding 
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g|“  ■  (x  sin  x)/(8in  x/x  -  cos  x) 


2 

X 

1 1 


^  + 
3!  - 


2  4  2  4 

1-^+2-  _i+2L..i:  + 

3!  5!  *•*  '  9\  At  * 


J-Jl 


2  -  2 
-X  +  3* 


=  5 


2!  4! 


which  is  va,lid  to  k^r  <  2»10"\  Thus,  we  obtain 


I  o 


A^.xp(-ikjr_,)(k^r2  -  5ikjr^  -  3)  =  Cexp(.ikjjr  )(k^jr^- 3^0- 5) 


2  2 


(3  +  5ikjr^  -  k^r^) 


(3  Jikjjr 


2  2v  ‘ 


Since  kjjr^<^1,  quadratic  terms  are  negligibly  small 

with  respect  to  unity  (at  an  accuracy  of  if).  Hence, 


0 


./  lO  N/,.'!!  I'O  N 

- - - )(1  +  - -  +  ...  )  = 


1,2  2  ,  ,  2 
k  r  -  k  k^-.r  +  ... 

-  A  (1+  -1^-2 - L-IL^ -  ) 

1  +  ikjjr^  +  ... 

Only  terms  of  the  second  order  yield  a  correcting  contribution. 
^OT  kjr^<tlO"  and  kix^o<10"  ,  i.e.  f  <100  kc/sec  and 
r  <300  m  this  correction  is  smaller  than  if-  and  thus  far  below 
the  reachable  accuracy  of  measurement. 


2 , 3  Maximum  moment  of  magnetic  antennas 

As  bhe  strength  of  the  excitation  of  electromagnetic  waves  in 
any  point  r  is  proportional  to  the  antenna  moment,  we  studied 
what  maximum  change  can  be  reached  by  filling  the  antenna  coil 
wiih  an  iron  core  ox  a  similar  material  of  high  permeability. 
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From  the  well-known  relation 


m 


V-  _ 

1  +  (u  -  1)N 


m 


o 


^oln 

where  N  is  the  demagnetizing  factor  and  m^  =  — j —  is  the 

moment  of  the  unfilled  antenna,  it  follows  by  a  limiting 
process 


m 

,  .  o 

lim  m  =  — 

p-is^OO  i'* 


i.e.  that  no  more  than  the  l/N-fold  of  m^  can  be  reached  by 
increasing  the  permeability.  As  m^  depends  on  I  as  well  as  on 
the  geometrical  dimensions  of  the  coil,  every  coil  filled  with 
material  can  be  replaced  by  an  air  antenna  with  the  same 
moment,  if  the  dimensions  and  number  of  turns  are  chosen 
adeauately. 


2,4  Sub-scale  tests 

The  possibility  of  sub-scale  tests  in  the  laboratory  was  taken 
into  consideration  as  the  theory  studied  so  far  requires  certain 
conditions  as  to  the  geometry  of  the  range  of  measurement  and  its 
physical  properties  (e.g,  unbounded  space,  homogeneity,  etc) 
which  are  difficult  to  fulfill  in  practice.  In  our  case, 
however,  such  experiments  proved  impossible  mainly  because  of 
the  required  high  conductivity  of  the  model  and  its  frequency 
dependence.  A  survey  of  the  references  in  this  connection  has 

tr£sv\  •iiy\  fPTiT  *1=  I O  i 

idWAA  A-At  I  |  .  9 
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2,5  Determination  of  rock  conductivity  from  VLF-propagation 
measurements 


In  the  propagation  of  VLF  waves  the  conductivity  of  the  earth 
plays  an  important  role  for  the  attenuation  of  these  waves. 

The  conductivity  was  determined  by  extensive  calculations  for 

’—A  —'I 

a  large  range  of  all  variables  (conductivity  o  =  10  ~  5*10*'  •”  m” 

frequency  v>  =  1  cps  -  100  kcps 

distance  r  =  100  -  20  000  ra) 

The  formulas  used  were  derived  with  the  aid  of  a  Hertz  potential 
for  a  dipole  in  the  dissipative  medium.  (Exact  derivation  see  [I'l) 


I  H  I  cB  f  ( r ,  \)  j  o )  = 


“  1 T  +  5 

r  r^ 


u  a 
0 


2"  ’^0°' 


1/2 


j.  v. 

^4  ( 


a  =  t,  •  o  =  2  tinJ  •  o 


For  purposes  of  clarity  some  of  the  values  calculated  are 
compiled  in  Fig.  2,1.  The  conductivity  of  the  medium  is  to  be 
concluded  from  the  value  of  the  attenuation.  In  order  to  separate 
two  values  of  0,  the  corresponding  value  of  B  must  have  a 
difference  of  at  least  5/  because  of  the  measurement 

inaccuracy.  Only  thus  approximately  safe  conductivity  values 
can  be  obtained.  Fig.  2,1  shows  that  only  at  r >  1000  m  a 
differentiation  of  the  field  strength  received  can  be  expected. 

At  a  frequency  of  1000  cps,  only  from  '^lO  000  m  upwards  a 
conductivity  up  to  '^'5*  10  can  be  determined  from  the  attenuation. 

-5 

At  higher  conductivities,  approximately  from  10  upwards, 
measurements  give  reasonable  results  already  at  5OOO  m.  At 
500  cps,  however,  and  a  =  10~^  a  measurement  is  expedient  only 
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at  a  distance  of  20000m,  i.e.  a  difference  of  is  observed 

-4  -5 

with  respect  to  neighboring  values  of  about  3*10  and  5* 10  . 

-7 

For  a  =  5*10  a  distance  of  20  km  is  necessary.  Table  2,2 
shows  the  minimum  measurement  distances  for  o  for  the  ordinary 
frequencies  on  the  assumption  that  a  change  of  the  oonductivit5^ 
by  a  factor  of  3  results  in  a  change  in  B  by  about  57  •  At  a 
given  measurement  inaccuracy  of  5'  this  will  probably  lead  to 
a  reasonable  determinability  of  the  d. 

Table  2,2 

frequency  minimum  distarce  for 


N) 

Q 

II 

0 

1 

a  =  10“"^ 

a  =  10“^ 

a=  5.10"’^ 

500  cps 

1500  m 

2000  m 

10  000  m 

20  km 

1000  cps 

1200  m 

1500  m 

5  000  m 

10  km 

3000  cps 

800  m 

1000  m 

3  000  m 

8  km 

lOkeps 

300  m 

V  . 

0 

0 

1  500  m 

5000  m 

100  keps 

100  ra 

150  m 

500  m 

1500  m 

The  measurements  have  hitherto  been  made  |1  ]  at  distances  belov/ 

1000  m  (in  exceptional  cases  up  to  15OO  m)  at  5OOO  cps.  Since 

-4 

the  conductivity  is  very  likely  belcw  10  no  definite 
conductivity  can  be  expected  from  these  measurements  according 
to  Table  2,2.  To  confirm  this  theory  various  measurement 
series  from  [l j  with  widely  differing  o  are  interpreted. 

In  order  to  compare  the  theoretical  values  (f(r,a))  v/ith  the 
experimental  results  (P)  the  mean  value  of  the  c(a)  is 
calculated  for  different  a  (different  a)  from 

C(a)*f(ra)  =  B(r) 

for  constant  a  and  the  corresponding  curves  are  compared. 
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This  comparison  is  limited,  however,  to  the  slope  of  the  curves. 
No  'absolute  comparison  was  made.  The  relative  deviation  of  the 
C(a)  from  the  mean  value  C(a)  is  the  measure  for  the  goodness 
of  the  approximation.  For  measurement  values  B  see  fl'j 

Table  2,3 

Measurement  Gertraudi 


a 

10-« 

O 

1 

>07 

5‘10"^ 

10“^ 

1 

O 

• 

GO 

8-10 

Ag/c 

0.7 

0.25 

1.2 

>  4 

>4 

0.7 

Measurement  Gertraudi 


a  10“^  10“^ 

5*io“5 

o 

1 

Ac/C  0.5  0.5 

0.5 

0.5 

Measurement  Lafatsch 

a  lO"*^  10“^ 

10“^ 

io“^ 

CD 

* 

o 

1 

i^c/C  0.6  0.8 

0.7 

0.6 

>5 

Measurement  Schwaz 

a  8*10”^*^  10”^ 

10"^ 

10"^ 

1 

O 

• 

GO 

Ac/C  0.8  1.1 

1.0 

>4 

This  clearly  indicates  that  the  conductivity  cannot  be  determined 

-4 

from  these  measurements.  High  conductivities  above  10  can  be 
excluded-  definitely  which  is  in  agreement  with  the  data  in 
Table  2,1 


Effsct  of  the  displacement  current j 


(for  exact  theory  see  Ml)  For  a  conducting  medium  the  wave¬ 


number  becomes  complex  -  ikg*  ''^ith  the  measurements 

hitherto  made  a  conductivity  of  ^  5*10  m  is  likely. 
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Table  2,4  gives  a  survey  on  the  error  arising  due  to  the 

*"  5 

neglection  of  the  displacement  current,  at  o  =  10 
Table  2,4 


frequency 

minimum 

distance 

difference  between 
k^  and  k^ ;k2 

difference  : 
B  and  f 

5  kcps 

2  km 

1 

'  t 

1^  kcps 

1 .5  km 

ro 

1 

5^  kcps 

1  km 

5 1- 

2.5f 

100  kcps 

500  m 

20  io 

10 

and  kg  give  the  exact  values  of  the  wave  number  with 
displacement  current,  for  k^  the  displacement  current  was 
neglected.  Since  also  the  measurement  distance  enters  the 
error,  the  error  becomes  larger  for  smaller  a  and  correspondingly 
smaller  for  higher  conductivities  because  the  measurement 
distances  required  become  larger  or  smaller. 

Since  the  measurement  inaccuracy  is  at  about  it  is 
justified  -  because  of  the  very  short  measurement  distances 
available  -  to  make  the  measurements  at  10  to  50  kcps  above 
1 300  m  without  taking  account  of  the  displacement  current  in 
order  to  be  able  to  determine  the  order  of  magnitude  of  the 


conductivity. 
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2,6  Horizontal  magnetic  dipole  at  the  depth  h  under  the  earth's  sur¬ 
face  -  Calculation  of  the  radiation  field  by  the  integral 

method 


1 )  The  radiation  field  of  electric  and  magnetic  dipoles  in 
vertical  and  horizontal  position  on  the  earth's  surface  and  an 
expansion  of  the  field  strengths  in  powers  of  the  numerical 
distance  was  given  by  Sommerfeld  [9]  s-  similar  manner, 

Ott  [11]  calculated  the  Hertzian  vector  if  of  a  vertical  electric 
dipole  located  in  a  non-conducting  but  optically  denser  medium. 

By  a  proper  path  of  integration  he  succeeded  in  separating  the 
excitation  into  a  main  wave  and  a  secondary  wave  called 
flanking  wave. 

In  the  present  paper,  Ott's  method  is  used  to  calculate  the 
radiation  field  of  a  horizontal  magnetic  dipole  within  conducting 
earth  under  a  plane  earth's  surface.  Ott's  separation  into  main 
and  secondary  waves  was  shown  to  be  feasible. 


2) 


The  vector 


air  j 
i 


and  the  boundary  conditions : 

A  horizontal  magnetic  dipole 
„  be  at  the  depth  h  under  the 
interface  earth-air.  The 


Let  the  conductivity  of  air  be  zero, 


electric  properties  of  the 
earth  are  characterized  by 
the  conductivity  a,  the  di¬ 
electric  constant  £,  the 
permeability  u  =  and  the 
wavenumber  kg . 
its  v/avenuraber  k  =  U'/c. 
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the  following  is  to  hold: 


k  =  kgH 


Por  a  time  dependence  expf -iotj 


2  2 

.2  «  .  4  ^ 

"  2  ^  ot  -i  ■'  2 

c  ”^o  n 


|n|<1 


a 

n  = 


2  o 
£  + 


2  2 
^o 

.-12.- 


(jj  angular  velocity,  c  light  velocity,  «  8.859*10"  [As/Vm^ 
kg  and  n  are  mostly  used  in  the  following  form: 


kg  =  Ikj,|  e 


1<P 


n  =  |n|  e‘^‘<' 


The  electric  and  magnetic  fields  are  obtained  from  a  two- 
component  Hertzian  vector  tT  =  (l?  ,  0,  H  )  (cf.  Sommerfeld  [9]). 
By  suppressing  a  optional  constant  and  the  time  dependence  here 
and  in  the  following: 

H  =  k  Tl  +  grad  div  T?'  z  <  0 


At  the  interface  z  --  0,  Tl  must  satisfy  the  following  four 
conditions : 

IT- 


2r 

i-Q  =  n  n 
Ex  X 


all 

Ex  X 


dz 

a7i„  an 

Ex  X 


dz 

an 


ax 


dx 


_ z 

dz 


air 


(1) 


Ez 


dz 


3)  The  x-component  of  the  excitation: 

n  is  composed  of  a  primary  excitation  IT  and  the  reflected  wave 
Hi'X.  o 

iJ  .  No  primary  excitation  is  assumed  to  exist  in  the  air.  Thus, 
S 

n  is  equal  to  the  refracted  v/ave  11 
X  s 

^Ex  “  ’’'o  ^  "r 


■n 


g 


e 


Ti_  =  —  can  be  represented  in  the  cylindrical  coordinates 

^  K 


r,  z,  <p  i-^s  a  superposition  of  pl-^ne  v/aves  in  integral  form; 


ik 


tf  a 
0 


E  I  j  ik_(r  sin-  cos((p'-(p)  +  (z+h)  cos  §) 


i  *‘'-E 
e 


2-^  ^  sin'')d  (2) 

In  the  exponent,  the  plus  sign  holds  for  z  >h,  the  minus  eigr> 
for  S\  h.  The  integration  with  respect  to  tp*  is  carried  out  from 
0  to  2n,  The  integration  with  respect  to  is  conducted  along 
the  curve  C'  in  the  complex  '^-plane  (Pig.  2, ’h).  Likewise,  -che 
reflected  and  the  lefracted  wave  can  be  represented  in 
integralf orm : 

ik-,  ff  ik^i'r  sin  .>  coe((p' -9)  +  (z-h.)  cos'<M 


f(>:^)sin  >d  Mip' 


ik_  ,  iv  !  r  sin '7^  cco(9'-9)  -  (z-h«)cos  i7pj 
ff„  =  ~i|e  ■  2  ^  g('>)sin-^dvd9’- 

The  amplitude  functions  f('')  and  g(i!?')  are  determined  from  the 
boundary  conditions  (1) 

m  .  -1 .  n2e{«  s(p  =  2 . 


C  /  a  .  i  f  N 

n  cos  X7+  (n  -  sin  v) 


■m 


2  /  2  .  2qs1/2 

n  COE  (n  -  sin  u )  ' 

2  . 

-  r.T  n  iV  )  ' 


n  C03W  +  (n  -  Gin  p ) 
Tor  this  v/e  find 


k^,h  cesv- 
£/ 


-kgh^  cos’^  =  kh^  cos^'^g 


and  the  law  of  refraction 

kgSini^  =  k  sin’*^ 

Integration  in  (2)  and  (5)  vrith  respect  to  9'  gives  the  Bessel 
function  J  (k„r  si.n'fi)  under  the  integral  sign.  Finally,  the 

O  Xi 

follov/ing  is  obtained  for  the  x-component  of  T  using  the  relation 

J  i  (eI  +  n^)  and  H^(-x)  -  hVx) 
o  20  o'  0  o'' 

ikj,(z+h)cos3- 


ik, 


rE 


- 4.  — 1  '  H^(k  r  cin>^)  e 

R  2  ■  E 

C 


f(-;>)8ini^d  <>  (4) 


n  -  15 


ik  .  -ik^.z(n^-sin^i‘^'.)^'^‘^+  ik„h 

T  _ h  Til /"i  ■  •  ^  Si  E 

,i^  =  -  .  H^(kgr  sin.')  e 


gO^)sin.i'd  '.(4) 


The  curve  C  is  shown  in  Fig.  2,b. 


^  I  s/ 

H  , 


•  '  9 

— _> - nr^-i  (.  v 

‘t  \  i  . 


r  A’p 


-Tt " 

V 

/ 

/ 


F  :C  ■:< 
*  ?  ^ 


-A’'  N  \ 


\  ./C.i 

'A  I  I 


- 4 


:/  ■  \  V  ■.: 
\  i;  ■ 

f  :  • 


4)  The  z-component  of  the  excitation: 


L  -  ^  Fig.  2,b 


Here,  also  has  no  primary  excitation.  The  boundary 
conditions  demand  a  slightly  modified  formulation  (cf.  Sommer- 
feld  [9']) 

ik—  n  f  .  ik-[j>(  z+h)cos  u 

^zE  T"  ^  ^o^V  si"'"')®  s(;5)sin<'dV' 


^  a  ^  oos'^ 

=  coscp  —  —  H^(kgrsin,/)e  s(,V)  gin  ?d  J 

^  C 


i(-.5)  = 


2  cos  >H(n^  -  1 ) 


k_[n^cos'^+  (n^  -  sin^i*?)^^^j  l_cos  i>  +  (n^  -  sin^u)^/^  ( 


5)  Calculation  of  the  integrals  by  the  saddle-point  method; 

2  2  C  1  /2  r 

The  root  v;  =  (n  -  sin  ^  '  is  contained  in  the  functions  f  ( '') 

and  s(t^).  To  obtain  w  as  a  utiique  function,  v/e  choose  the  real 
axis  of  the  w-plane  as  branch  cut  and  image  it  in  the  .^-plane. 
Branching  points  appear  in  the  periodic  interval  [-irjnl  at  sin^'^=+n 
and  sin^Tt-.*^)  =  +n  (Fig.  2,b).  The  equation  for  the  branch 


II  -  14 


cut  Iin(v/)  =  0  reads 

sin2^3'^  sinh  2^  =  -2ln^!sin2<p  =  const. 

It  describes  a  curve  through  the  branch  point;  its  solid  part 
in  Pig.  2,b  renders  the  axis  Im(u)  =  0  and  its  dash-lined  part 
the  axis  Re(w)  =0. 

f('})  and  g(\}')  have  poles  for  cos  =  -w  /n^,  \7ith  the  solution 

P  P 

sin  =  +(n^/l  + 

2  2  1  /2 

In  the  pole,  the  root  has  the  forn  v/^  =  +n  /( 1  +  n  )  ”  Y/ith 

T  .  2  1-  r  2,  -i2p.  ,  ,  2//  2xl/2  K-i(2ip-p) 

1+n  =1+n  e  vr  reads  w  =  +  n/(1+n)'  e 

p  p  — '  '  '' 

vrhere  2«p-p>0.  This  shows  that  in  the  upper  Riemann  sheet 
(lni(w)>0)  the  negative  sign  Jir.st  be  chosen  for  v;^.  Thus, 


cos 


S'  =  cos  coshif2  ~  i  sin’5'^  sinhi^2  =  !''/('• 


This  gives 


;  •;^2 


2-  >2^" 


P„: 


0 


in  the  upper  sheet.  In  the  lov^er  sheet  (im(vir)  <0),  the 
follov/ing  is  found  analogously ; 

P,: 


■\>o 


P4:  1>2>0 


-  I  <'■  <  Tt 

-  n  4  f 


Besides  it  can  be  shown  that  the  poles  alv/ays  lie  between  the  real 
axis  and  the  curves  given  by  .sin  2w  sinh  2^„  =  const.  (Pig.  2,b) 


In  the  0  point  of  the  upper  sheet  has  the  form  ;|;n= +|n|o 

As  in  the  upper  sheet  Ira(vf)>  0,  the  negative  sign  is  applied. 

Thus,  Wp  has  a  positive  real  and  imaginary  component  before  the 

branch  cut,  the  signs  only  changing  beyond  the  cut. 

2 

In  we  replace  by  -1  +  n  The  contribution  of  the 

i^X 


-icp 


summand  -1  can  be  immediately  found,  and  we  obtain 

+  I, 


ik,.,R 


ikgh  ‘ 


Ex 


R 


'1 


11-15 


ik^n"  !  ^  ikg(z+h)cos''^ 


(k„r  sin  %)e 
2  E  ' 


g(’5)sin  ’"^d 


W 


To  calculate  we  lead  <3,  the  path  of  the  integration,  somewhat 
out  of  zero  and  first  replace  H  (k„rsint>)  by  its  assymptotic  vcluo 

O  ill 


„1/,  .  2  ^  E 

H  (k_,r  sin  j  =  if  “1 - : — ^ 

o'  E  ■'  )  uk^rsinA 


.  .  Q  7t  \ 

>  i;,k„rsini^‘ -  j  ) 


The  eligible  exponent  is  now 

i[kgr  sin’i  +  kg(z+h)  cosi^j 

or,  after  introducing  polar  coordinates  a,  R’  in  the  r,z-plane 


r  =  R'sina  z+h  =  R'cosa 

i[kgR*  sina  sin^  +  k^R'  cosa  cos  =  ik^R'  cos(a-i*i) 

The  real  component  of  the  exponent  becomes  <0  for 

-K  +  -^  +  a<  ^  +  a  ‘^2  >  0 

a<6^<^~+a  ^2^® 

The  path  of  integration  may  run  in  the  hatched  part  of  Pig.  2,b. 
Since  ihgR'  sin  (a--^)  =  0  for  a  =  i^,  the  saddle  point  lies  on 
the  real  axis.  The  falling  line  Im(cos(a-i5) )  =  const,  runs  near 
the  real  axis  from  the  pass  a  under  155°  toward  +(;^  -  ioo)  +  a. 
At  a  greater  distance  from  the  pass  its  gradient  deviates  from 
the  value  -1  but  the  pass  road  is  no  longer  forcing  there.  V/hen 
we  transfer  the  path  of  integration  into  the  falling  line,  the 
branch  cut  disturbs  us  in  the  second  quadrant  for  a>a^=arcsin  n. 
Prom  the  end  point  ^  +  ioo  no  direct  path  leads  to  the  pass. 
Therefore,  we  must  go  round  the  branch  cut  and  take  the  foJLlovnng 
path  for  this  purpose;  from  —  +  ico  on  the  upper  sheet  to  +ioo, 
from  there  along  the  branch  cut  to  its  origin  and  on  the 
other  side  back  to  +ioo,  from  there  on  the  lower  sheet  to  the 
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end  point  of  the  pass  road  -  +  a  +  ioo  and  from  there  over  tne 

pass  =  a  to  the  old  end  point  ^  -  ioo  (Fig.  2,b,  curve  6^). 

As  the  integrand  on  the  approaching  paths  vanishes  at  infinity 
the  integral  I^  is  split  up  into  two  parts  Ij^  and  Ilj.  I|  re¬ 
presents  the  saddle  value  I||  the  integral  along  the  branch  cut. 

In  Ij  we  replace  g(d')  by  the  saddle  value  g(a),  which  is 
permitted  for  a  and  draw  this  factor  out  before  the  integral. 

If  we  use  the  original  value  instead  of  the  asymptotic  one  for 
the  Hankel  function  under  the  integral  sign,  we  obtain: 


ikgn^g(a)  ^ 


.  ik„(2+h)cos  1*^ 

Ho(kErsin\^)  e  sini?dv^  = 


=  n  g(a) 


k 


or,  according  to  (4) 

Analogously  we  obtain 


ilCpR 

^  +  f(a) 


ik^R' 

y 

_ 


+  I" 


TT. 


=  cos(p  I-  [s(a) 


+ 15] 

!> 


Ez  ^  c>r  R» 

with  a  corresponding  definition  of  I^. 

Near  a  =  •^,  where  g(a)  and  s(a)  disappear,  the  approximation 
becomes  useless,  'Ve  must  develop  g(ii)  and  s(i5)  in  1?  =  a  and 
take  account  of  further  terras.  Confining  ourselves  to  the 
second  order  term  we  obtain 

2 
J 

^E" 


T,  r  2  /  \  in  g"(a)-|  e _ 

I^  =  [n  g(a)  -  g,  E,  J  El 


,  “i-O  n  ^ 

i2i 


u 


^2  ar 2k„R'  ■'  R'  J 

This  is  the  separation  ij.to  main  and  secondary  waves  as  shown  by 
Ott  [ll'l,  I^'  and  Ig  are  the  corresponding  quantities  in  our 
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express ions.  Marking  the  hitherto  calqulated  integrals  bj  th' 
superscript  M  ire  obtain  for  the  main  ware  in  the  earth  t 


n 


Ex 

Si 

E2 


ik^R 

V-  * 


? 


cos<p 


ar 


[s(a) 


ik^E* 

L _ _  1 

Rt 


n 


H 

Ex 


e 


ifcgR 


ikgR* 


.  2 


2fc^R- 
ikgR» 


n 


St  -  -  --yi)  -gT-  J 


E 


€)  fhe  secoffiiay  wave  I?  f-ivI  I^t 

_,  -  I  /. 

As  on  tlfe  hraneh  cut  (n  -  nin^'S) 


r/i 


+  't  »  where  t  is  a  real 


naaher*  we  bring  1^  with  the  aid  of  this  substitution  equation 
into  the  fora 


2  2 
i-nSr 


»G0 


2t 

t^n 


and  for  the  Integration  we  again  pass  over  to  the  complex 

2  1 

t**plaae#  Branch  pointo  lie  at  i’  -  +n  and  «  +t(l-n  )  if 

2  2  1  /z 

we  again  choose  the  real  aris  of  -  {n  - t  )  '  ,  the  equattor 


for  the  branch  cut  reads 


Tg  =  1 


i  r  ~  I  s^2n 
- 


e, 


?he  poies  are  given  by  +Iu^/(1  i  , 

they  lie  In  the  second  and  fcrtL  quadrants  of  the  t -plane  and 
do  not  disturb  us  when  we  tTans^ei’  the  path  of  integration 
into  the  pass  road. 

With  use  of  the  asymptotic  reprecentatron  for  the  llankel 
function,  the  exponent  in  questian  becomes 


^ V '  [:  sina(n^  -  +  coaa(  1  -n^+  ■?f  « 

^  .  0  i  ^ 

(  5  .  -  alh2a)  V2 

Thus,  the  saddle  point  in  question  lies  at  T  «  0.  The 
expansion  of  f(T)  in  r  »  0  reads* 

f(T)  ■  ikgR'cosCa-n'^^)  -  jqle^^T^  + 

ik^E'ein(a-ii) 

q  ■  — •-■  ‘  "  '■  * 

2Din»5Q  cost^Q 

The  equation  for  the  paasroad  ic 

'  *^2  *  ±Cc^  +  1)  t-|  c  -  cte0 

and  its  gradient  near  T  *  G  is  given  by 

tg-»]«*-tg^  0*  argfqj 

1/R(t)  •  i/fl'  +  n^(l  -  n^  +  expanded  in  T  ■  0, 

i/n(t)  =  i/h(o)  -  t/jj(o)^ 


f(r) 


•  •  a 


n 


gives 


where  only  the  second  term  makes  a  contribution  to  Thus, 

we  obtain 


ly  - 


■t 


lkj,R'co»(a-4H|  !  . 

mm  d  ' 


nhr-n^) 

and  finally 


e 


J 

-£ 


it: 


s  2e 
Ex 


iY 


T^+g-i^^otga 


3/2  I  -n^  J 


w  _  3  b 

Y  4  +  3<P  - 

and  analogously 

^  ik;p^nr+  (z+h:)il  'I 
f ''  -e_;^ _ 

!  r^ 


ri 


Ez 


cos<p 


.4rAi 

~  n  ± 


u2//  2x 

k  (1-n  ) 


o 

5r 


j*‘ 


l  ctga 
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^Dhus,  the  vector  n  is  determined  in  the  earth.  ‘Te  should  Tike 
to  mention,  that,  in  agreement  with  Ott  [11 j,  the  secondary  wave 
disappears  in  the  angular  space  ot  <  =  arcsin  n  since  the  branch 

cut  in  the  fourth  quadrant  must  he  also  avoided  and  thereby  I!j* 
and  I^  ere  Just  eliminated.  Besides,  this  representation  is 
identical  to  that  originalTy  given  by  Sommerfeld  [9]  [lOj;  if 
kj,sin»9in  (j)  is  replaced  by  X  the  expressions  go  over  into  the 
Sommerfeld  integrals.  Integration,  however,  is  different  from 
that  in  Sommerfeld 'e  paper  and  docs  not  give  rise  to  any  residue 
waves  since  the  path  of  integraiion  does  not  include  the  poles. 


7)  gho  vector  -g  f jr  g  <  0 ; 

By  means  of  the  transformation  sini^  *  n  sin^  we  transform  IT^  into 

-ikz  cos  j  +  ikgh^l  -n^  sin^| 


^  ffi^Ckr  8  inf  ) 
p(f)  * 


p(f Jsinfdj 


2n  cost 


21^  2  .  2c'  i 

n  \|  i-n  sin  F  +  ncos^ 

Calculated  analogously  to  the  former  integrals  and  using  the 
polar  coordinates 

r  =  Rq  sinp  -  z  *  cosp 

this  gives  for  the  antenna  depth  h  =  0 


ikR. 


11 !!  =  (p(p) 


ikR 

V-S _ - 

zkn  f  r  ; 


o  o 
ikR 


-n  a  ,  .  .  s 

•*2  =  C0S(?-^1S(- 


^  COSV^  ^ 


i  3”  (  P  )  ^  £. 


ikR 


o  -1 


2kR_  '  R 


P  0 
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and  for  the  secondary  wave  5 


nS  =  J!e1_( - 2 - 

*  2(l-n  )  \  k_R  elnpiq-r/ 


\  i''siSoL=^"P  oosp\|n^-l]+  lYj 

5,1  '  (5) 


COS(p 


2  ^ 
n  +1  a 


\  ® i»ip-cos 3\ n^-1  J+iCYg-  f )' 


2ln2-1  k_R  sinp  Iqpl  ^ 


Y  -  I  .  i£ 
^242 


ik„R 
•h  o  r 

““T“  [sinp  - 


r  =  arg[q2l 


The  formulas  (8)  and  (9)  give  n  for  the  case  h  »  0.  For  finite  h, 
this  representation  may  he  taken  as  a  good  approximation  for  the 


following  case: 


or  on  the  pass 


kgh^  1-n^ain^^  kR^co8(j.“p) 


1  -  n^sin^p 


(10) 


As  soon  as  this  relation  does  not  hold  any  longer,  the  above 
integrals  can  be  calculated  in  the  same  manner  for  finite  h,  while 
the  saddle  point  is  shifted  by  a  factor  depending  on  h.  Concerning 
the  range  of  validity  for  the  representation  by  the  saddle  point 


method  we  should  compare  the  theoretical  study  in  [l 2],  according 

to  which  the  neglect  in  an  integral  of  the  form  dz 

is  of  the  order  i/m,  which  means  in  our  case  that  the  error  becomes 

of  the  order  |  l/kj,R'|  ,  The  expression  fo:;  n  can  be  ciraplified  if 

we  leave  the  hitherto  observed  generality  and  restrict  ourselves 

to  n^lO”  ,  In  the  exponent  of  (7)  the  term  (:1-n  sin  '  can 

ikgh 


be  considered  constant  and  the  factor  e 


can  be  drawn  out 


before  the  integral .  That  means  bhat  we  have  found  a  representa¬ 
tion  also  for  large  h,  for  which  relation  (IO)  does  not  hold. 

The  expression  for  h  =  0  has  only  to  be  multiplied  by  the 


ikgh 


,  Thus  the  effect  of  h  is  characterized  by  a  phase 


factor  e 
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shift  and  a  damping. 
Hence, 


ik„h 

e  P(p) 


-ikgh 


=  e  cos^p [s(p) 


and,  analogously,  the  other  formulas. 

8)  The  field  strengths; 

Although  the  calculation  of  field  strengths  in  the  earth  is  a 
little  complicated  it  does  not  involve  very  great  difficulties 
for  numerical  computation.  Since  in  (6)  the  first  term  contains 
R  while  the  other  terms  are  written  in  R',  we  transfer  the  center 
of  the  coordinate  system  into  the  mirror  point  and  calculate  in 
the  spherical  polar  coordinates  R',  a,  <p. 

Because  R  =  R'  -  ^R'/dz  =  R'  -  2h  cosa 

i/r  =  i/r  -  .^i/R‘/^’z  =  1/R'(1  +  2h  cosa/R') 
we  can  write 

ik„R  ikpR'  -i2k„hcosa 

e  /R  =  e  /R'(1  +  2hcosa/R')e 

The  functions  f(a)  and  s(a)  can  be  simplified  for  a  >20°  and 

nCIO"^  to  f(a)  =  -1 

s(a)  =  “2/kg  ctga(cosa-i  sina)  =  s^(a)/kg 

Thus,  we  obtain  the  components  of  the  main  wave 

after  long  intermediate  calculations  ir.  the  following  form. 


where  A  is  an  antenna  constant; 


=  A  G03M>  e 

n  * 


i  (kgB '  -  u.  t )  ^  g^  ( a,  h ,  kg  ) 


1  ’  ' 
iii  R' 


=  A  cos<p  e 


fc(kgR'-k.t:)  ^  f.(a,h,kj,) 


(11) 


=-A  sincp  e 

<p 


i  (k„R  •  -  ...it)  A  h .  (  a ,  h ,  kg  ) 

z  - 5 — 
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2 

*  kg  (i  sina  cosa  8^(a)  -  sina) 
gg  •  (-2hkgCOsa  sina  -  ikg  sina  e 
^  »  ikg  Rcosa 

gj  =  sina  e"^(-14h  cosaikg  +  2  -  4h^kg  cos^a)  + 

+  2is^ (a)  cosa  (—r —  +  sina  +  i  cosa) 

•  sina  ' 

^4  *  [sina  e"'^(-24ikgh^  cos^a  +  18h  cosa)] 
g^  -  [sina  e“^(48h^  cog^a  -  8h^  cos^aikg)] 
gg  =  sina  e  40h^  cos'^a 


^  sin^s^(a)  +  cosa(e“’^-1  )J 


^2  ^c^  ^  kgS^(a)sin^aJ  + 


*  ^  *  2  ^  2  * 

+  e  ^[2hkg  (cos  a-sin  a)  +  ikg  cosa] 

^5  “  I®  ^(4Rikg  cos^a  -  cosa  -  Sk^h^  cosa  sin^a  -  6ikgh  sin^a) 

*"  cosa  (sina  +  -r^ —  +  i  cosa) 

t’a  1  sina  ' 

^4  =  e  ^  (4ikgh^  cos^a  -  6h  cos^a  -  8kgh^  cos^a  sin^a  + 

2  2  2 
+  6h  sin  a  -  20ikgh  cosa  sin  a) 

-1/22  2  32  2  ox 

=  e  (24h  cos  a  sin  a  -  24ikgh^  cos  a  sin  a  -  12h  cos^a) 


(e-  ^  -  1) 


_  £  2 

hg  =  e  ^  (2kg  h  cosa  +  ikg) 

^5  “  6  ^  (4h  cosa  ikg  -  l)  -  s^(:a)  +  sina  +  icosa) 


^2  2 

e  ^  (4h  ikg  cos  a  -  6h  cosa) 


h.  =  e” ^  (-12h^  cos^a) 


e  ^  (-8h^  cos^a) 
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The  formulas  for  the  earth’s  surface  z  =  0  are  much  simpler. 

In  a  polar  coordinate  system  r,\}',  (p  the  origin  of  which  is 

on  the  surface  perpendicular  over  the  dipole,  the  following 
2 

holds  for  n  1 


=  A  cos»p 


=  A  cos<p 


=  A  sintp 


,.,2 

ikr 

C 

(  - 

lik  6 

1,5 

2 

r  2 

k 

r 

r 

5 

2ik.<e 

•b 

ikr 

e 

(  - 

lik  ^  6 

k^ 

2 

r 

r 

2i4e  " 

ikr 

e 

/,,2 

ik  2 

(12) 


r  2 
r 


where  the  dependence  on  time  has  to  be  included  yet. 

The  results  for  the  secondary  wave  in  the  earth  are  as  follows; 


.  ik  R’cos(a-6  ) 

2Ae  ^  e  f,  2/  .  „  2  .  x 

-  cos<p - ^ - 2 -  k  (,sina  -  n  cosa  -  2n  sinofi 

^  R'  sin'^a 


2ikg(cosa  +  4n  sina)  ^2 


R' sina 


2  2 
R'  sin  a  - 


1 


ik„R’cos(a-T}'  ) 

2Ae^^  e  ®  °  \ 

=  costp  — ^ - - - ^ — -  1  kglcosa  +  n  sina)  + 


a 


2  2 
R'  sin  a 


2ik_,  sina 

^  -c!  ^  12  cosa 


2  2 

R'  sina  R'  sin  a  - 


ik_R '  cos  (  a-1%)  „ ,  , 

2Ae  -  •  e  '1,2  - E"  4 

H  =  -Sintp  — r - 5 - x -  k^  +  - - — —  -  — 

^  R’'^  sin  a  '  ^  -.2  ^...2. 


R'sina  R'  sin  a 


sm 


=  n 
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and  on  the  earth's  surface  z  =  0 


»  C08(p 


cos(p 


H  «  -sin<p 
9 


ikgh+iY 

2iAe  k^ *  ikr 

E  e 

-ik^(l-2n) 

«» 

.kn!. 

12  in' 

k^ 

r 

2 

r'^ 

ik„h+iY  , 

2iAe  ^  ikr 

£i  8 

-ik^n(2+n) 

+  6>S!!' 

k4 

r 

ik  h+iY 

2iAe  ®  k;  ikr 

E  e 

Tt  0 

2" 

.  +  ilP- 

'  rs 

1 

(14) 


Comparing  (14)  to  (12)  one  can  see  that  the  of feet  of  the 

secondary  wave  increases  with  increasing  distance;  from  a  certain 

•  .  .  2  9 

minimum  distance  on  —  for  k  ft^lO  about  8«10  m  —  it  will  even 


exceed  the  main  wave.  As  the  phase  shift  of  the  secondary  wave 
Y  *  7  +  and  p  = - *  ® 


arg 


ikpR'3in(a-'fi'^) 


2  sin^  cos'^ 


follows  for 


“o  ’'o 

kg  «  A  (i  +  1),  (displacement  current  neglected)  for  n<  10”\ 
h  -  250  m  and  r>  2.5.10^  m  p  =  arg  [(l  -  i)n].  Thus  y  is 
independent  on  r.  V/ith  these  restrictions  the  amplitude  of 
the  field  strengths  at  z  »  0  is  given  by: 


H  I 

r*' 

A 


4A  -Ah  1 
cos<p  ^  e  -r- 


/r- 


Isl  ,  2  k^  1  6 

^  '  2a2  ^  r^ 


4k  6k^  1"*^ 
A  r  ”  ^2  ^2 


1/? 


H 


*^1 


cos(p  e 


Ah  1 


id  .  liSl  ± 

"  2 
c 

“1  2 


A  A^r  r^ 


(15) 


k^  .  k^  -  6k A^  1 


2A 


2 


^  id  e-Ah 

A  k^  r^ 


r 

/ 

k^  2 

2 

2 

k  2k'^ 

'1 

1 

.2  2 
At  r 

r  “  .  2  2 
A  r 

} 

2,7  Calculation  of  the  dipole  field  using  the  VLF  approximation 

The  range  of  validity  for  field  strength  formulas  based  on  the 
saddle  point  method  is  somewhat  restricted  by  the  fact  that 
some  parts  of  the  calculating  procedure  rest  on  the  assumption 
of  high  dimensionless  values  characterizing  the  distant 

field  of  the  dipole. 

This  makes  it  desirable  to  employe  separate  formulas  for  the  range 
between  the  remoter  parts  of  the  near  field  and  a  few  wave¬ 
lengths.  Sixch  a  separation  would  make  it  possible  on  the  one 
hand  to  calculate  the  strength  of  the  field  at  every  distance 
and  on  the  other  hand  to  compare  the  ranges  of  transition  as 
between  the  near  field  and  the  distant  field,  a  comparison  which 
is  of  particular  interest  experimentally  because  the  range  of 
the  measuring  transmitter  is  limited. 

The  starting  point  for  these  considerations  is  once  again  the 

electrical  vector  potential  n  which  is  related  to  the  field 

strength  IT  in  the  foliov^ing  way  ; 

H  =  k  7t  +  grad  div  n 
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Asaviming  the  same  orientation  of  the  dipole  as  in  (”>7) 

<'r‘  in  J5]f  it  has  two  components  ti  and  n  which  can,  in  the 
usual  way,  he  repi'^sented  as  an  integral  through  the  eigen¬ 


function  Jj^(X^)«exp  +  pa  *008  n^ 

00  00  _  ^ 

Earth  space:  «  a^J  ^  dX+Jfg(x)e  ^  J^(X^)dX 

o  0 

CD 

’'zE  “  J  J^(X^)dX»cos^ 

o 


Air  Space: 


\o  *  / 

o 

7  p  z 

^zo  “  i  ^0^^^®  °  J^(X^)dX»c089 


It  is  more  expedient  to  retain  the  Sommerfeld  formulation, 
integrating  in  the  X  field.  Here 


9  »  tan"  ^  x/y  ,  a^  =  I  4t^ 

*♦  <♦ 

The  conditions  for  the  continuity  of  E  and  H  in  the  boundary 
prane  a  ■  0  enable  the  amplitude  functions  f  and  g  to  be 


determined  as 
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,  .  ,  X  -'■e*'  "e'^o  -  >*04 

fpU;  *  -^  *  IT"  ®  — 2 - 7 

®  i^E  +  ugic: 


f  (X) 
o'  ’ 


=  A  •  X  •  e 


-Ugh 


2k 


E 


'‘o'^E  ^  '‘e‘'o 


(3) 


go(X)  -  ggC?^)  =  -  A 


V,  u2  ,2 

-Ue^  k^  -  kg 

X  e  - = - 


2 


E 


If  the  integrands  in  (2)  are  simplified  by  using  what  is 

called  the  VLP  approximation,  the  four  integrals  in  (2)  can 

he  calculated  exactly.  By  contrast  with  this,  the  me-thod  of 

calculation  explained  in  (2,6)  or  in  [5]  although  not 

neglecting  anything  contained  in  the  integrands  does  not 

lead  to  closed  solutions  of  the  integrals.  The  approximation 

2  2 

consists  in  neglecting  k^  as  against  kg  in  the  expressions  (j), 

which  is  permissible  as  long  as  the  condition  is 

weo 

satisfied,  the  latter  being  almost  identical  with  the  condition 
^  2 

(<— - )  $>1  which  has  to  be  satisfied  if  the  displacement 

"I  ^^5 

current  is  neglected  (and  which  is  valid  if  o/to  ^10  )  . 

The  vector  potential  Kg  in  the  earth  space  which  is  of  sole 
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interest  here  (the  horizontal  components  of  the  field  strength 
on  the  earth’s  surface  being  calculable  from  tL)  then  works 
out  as  fol<.ow8  using  the  relationship  ^  cos<p  •  : 

00 


xE 


0 


-Pg(z+h) 

e 


It 


-a 


zE 


B 


(4) 


-Pg(z+h) 


J  (\^)d\ 


The  physical  situation  postulated  by  this  mathematical 
approximation  is  to  be  interpreted  simply  in  the  sense  that 
the  equations  (4)  apply  only  to  such  fields  as  are  originated 
by  conducted  currents,  which  cannot  of  course  exist  in  the 
air  space.  The  theory  erected  on  the  basis  of  the  VLF 
approximation  is,  therefore,  a  direct  corollary  from  the 
calculations  made  in  the  second  report  ],  wherein  the 
displacement  current  was  suppressed.  Neglecting  the  dispiace- 


ment  currents  immediatclj 


O  -W  Wr-W.O 


the  range  over  which  the 


equations  are  valid,  limiting  their  application  to  the  distance 


within  which  the  conducted  currents  in  the  earth  space  are 
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due  mainly  to  the  primary  excitation  and  not  to  the  inductive 
effect  of  the  displacement  currents  in  the  air  space,  as 
would  he  the  case  if  k^R  were  large,  i.e.  in  the  distant  vacimm 
field.  The  resulting  expressions  for  the  field  strength  can, 
therefore,  he  applied  in  the  near  field  up  to  distances 
corresponding  to  a  few  vacuum  wavelengths  (X^  =  100  km  at 
f  =  3  kHa  [5000  cpa'l).  They  serve  to  complement  the  formulas 
of  the  saddle  point  method  in  the  short  distance  range. 

The  integrals  (4)  are  referred  hack,  hy  double  integration 
with  regard  to  a,  to  the  known  main  integrals 

CD  00 

e  ®  J  (X<?)  —  dX  or  ,  e  ®  J  (Xe)  dX 
o  o 

which,  hy  elementai’y  calculation,  work  out  as 


-Ik  Jb-  + 


2  2 

E  l  - 


or  ■  I^  (a)K^(  p) 


v;here  B  is  any  function  of  a,  I^  is  &  modified:  Bessel  function, 


ik_ 


is  a  modified  Hankel  function,  a  =  — —  (\j  B'^  +  -  B), 


ik^ 


^  I  —  — 

p  =  — —  (V  B  +  ^  +  B),  If,  as  on  page  Ii-10,  the 

2 


distances 


n  -  30 


+  {z-hf  and  R»^  =  +  (z+h)^ 

between  the  dipole  and  its  mirror  image  are  introduced,  the 
following  expressions  are  obtained  for 

E 

“  *E  \  R 

^z  "  ^ 


-ik^R» 

B 


R* 


V“K(p)^ 


and  the  following  for  the  horizontal  components  H  and  H,  of 

3C  y 

the  magnetic  strength: 


The  third  derivation  of  I^;(ar)K^0p)  for  z,  which  is  necessary 
for  calculating  and  Hy,  is  arrived  at  by  using  the 


recursion  formulas 


II  -  31 


dl  (a) 

=  I, (a) 
da 


4K„(0) 


=  -Kl(p) 


dl,{a) 


-  .  I  -  — 

da  ^  a 


<lK,(p) 


-K„(P)  - 


Ki(0) 


and  appears  as  follows : 


I  (a)K  (e)  =■  —  i  I  K  f  ik„  (1  - 

°  °  ^  2  j  0  o  t  E  R,  ^ 


4a 


(q  +  1)^J  +  I,K^  [kg  (1  -  q^)  -  ^  (q  -  1)5  . 


-  2  ^“k  -  ■')  (1  -  4  (1  -  4^) 


IK.  I  Ik  (q  +  1)  (1  -  q^)  +  k!  {q5  +  i)  + 
°  I  ^  aB  »  * 


(  q  +  1)5  +  -5^  q  (1  -  q^)  I  +  [ikg  ^  (q^  -  l)(3q-) 


j-  ^  (  q  -  1 )  (  q^ 

4  o 


1  )  -  i  ^  (q  *  1  )=  (q 

4  3 


1)1  ! 


On  inserting  (6)  into  (5)  and  making  a  brief  intermediate 


calculation  the  following  simple  expression  is  found  for  the 
horizontal  component  of  ^  in  the  direction  of  the  main 
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radiation  from  the  magnetic  dipole: 


-iM 


e 


1 


«xlx=o  =  ^ 


-iN 


1 


R'^  2 


“  j^-N^  +  iN^  -  +  6iN^  +  6-  2q^(l5+15iN^-6N^-iF^)j 


^  [l,K„  (51  -  15  q’  ^  q¥^) 


I  K  iN  (-1  +  3q  -  -  6q5  -  (1  -  q)(3-  ^5Q^)) 

I.K  iN  (-1  -  3q.q2  +  6q5- (l  -  q)(3  -  15q^) ) 

n: 

(-13q  +  •I5q^  -  q^N^)] 

with  =*  ^  ^)/E' 

Apart  from  the  VLF  approximation  nothing  else  is  left  out  of 
account  in  this  equation  for  H  .  A  glance  at  the  expressions 
included  in  brackets  within  trhe  second  and  third  summation  groups 
shows  that  for  h  — >  oo  and  for  a  finite  R  these  expressions 
tend  to  zero  by  reason  of  the  e-powor  vanishing  and  owing  to  the 
Hankel  function,  so  that  (?)  transforms  itself  exactly  into  the 
equation  for  the  field  strength  of  a  magnetic  dipole  in  a 
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homogeneously  dissipative  medium.  This  result,  as  well  as  the 


fact  that  the  field  strength  decreases  as  the  conductance 


increases,  is  trivial  as  it  can  at  once  be  inferred  from  the 


initial  equations.  It  is  possible,  therefore,  to  write  (7) 


in  the  form:  H  =  H  ,  +  H  , 

X  X  horn  X  suppl 


v/here  H  _  represents  the  additional  field  attributable  to 

X  suppl 

the  influence  of  the  earth-air  interface  behaving  symmetrically 


with  regard  to  2  and  h:  that  is  to  say,  independently  of  whether 


for  instance  the  transmitter  is  500  m  and  the  receiver  50  m 


belov;  the  earth's  surface  or  conversely  the  receiver  5OO  ra  and 


the  transmitter  50  m  below. 


Obviously  no  such  symmetry  applies  to  H  .  Simplifications 

X  ges 

can  be  made  in  (7)  if  the  distance  of  measurement  notably 


exceeds  the  height  covered  by  the  transmitter  and  receiver 


(q  — ?  0)  or  if,  on  the  contrary,  the  height  covered  notably 


exceeds  the  distance  of  measurement  (q  — *"1)'  These  two  limiting 


:a3S3  are,  hov/ever,  only  of  theoretical  interest,  since  only 


the  juterraediate  range  is  amenable  to  experiment. 


Tor  practical  evaluation  of  the  f'rrauJa  (7)  it  is  best  to  divide 
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into  a  real  and  an  imaginary  part,  as  the  apparatus  available 
at  present  allows  only  of  measuring  the  amount  and  not  the 
phase  angle.  Since  "the  arguments  of  the  cylindrical  functions 
arising  make  an  angle  of  155°  in  the  complex  planes,  I  and  K 


can  be  expressed  by  the  Kelvin  functions  her,  ber^ ,  bei,  bei^, 
her,  her^ ,  hei,  hei^  whose  argument  is  purely  real;  the 
functions  ber,  bei  have  the  argument  and  the 

functions  her,  hei  th"  argument is  (l+q.)E * k^/\i 2  ,  whore  has 
the  value  k,  =V — r-.  i!  short  interiaediate  calculation  then  leads 

j  H  d 

from  (7)  to  the  real  and  imaginary  parts  of  E  as  follows; 

‘A 


R  H  , 
e  xlx=o 


e”**  (Bp*cos  M  -  Bq  sin  M)  + 


+  -e''^(C.  cos  N  +  Cp  sin  N)  +  B  A.  +  C,(A„-Aq)  + 

,5  '1  2  '  0  1  3'  2 


+  ^5-^7  ^CAjCB^+Bg)  +  A^(B2-B^)  +  A^(Bj+B^)  + 


+  B^Ag  +  N(A^(B^  -Bg)  +A^^:CB^  ^ 


n  -  55 


M  = 

N  = 

q  = 

-^■1  =  ~  ("ber  hei  +bei  her  ) 

1  o  o  0  O' 

=  -  •?  (bei_^hei  -ber  her  ) 
i  o  0  0  o' 

^4  =  -  I  (ber^hei^+bei^her^) 
2^=3^“  i5q^ 

=  -1-q^-(6q^-5q) 

^2  =  (i+q)(3-i5q^)/2N^ 

=  -1-q^+(6q^-3q) 

=  -  N  +  1  +  I  -  2q^N  -  12q^ 

Cg  =  -2n2  -  H  -  (3/N)  -  (3/^2) 

C3  =  2nV 


k-R 

5 

k  R' 

> 

2+h 

R' 

A5  =  -  I  (bei^hei^-ber^her^) 
Ag  =  I  (bei^her^+ber^hei^) 
A?  =  -  f  (bei^her^+ber^hei^ ) 
Ag  =  I  (ber^her^-bei^hei^) 

B4  =  (l-q)(.3-15q^)/2N^ 

B5  =  -13q  +  15q^ 

Bg  = 

B^  =  2M^  +  M 

-  15qVN 

-  2q^N  +  (l5qVl^)  +  (l5qV^T^) 


The  result  of  the  numerical  calculation  for 
and  3n  can  be  seen  from  Table  2,5, 


J5iE  ta 

TO  ’  10 


f/i  a 


II  -  36 


Table  2,3 


R  = 

10 

lofTo 

10^ 

j'ra"! 

m 

1 

0 

3.i6*io‘^ 

1 .025*10"^ 

7.96' 10"® _ 

2.5202»10"^ 

i 

-1 

;iO=57l*10 

m 

10"^ 

3.1-10“^ 

1.009*10"^ 

7.78*10"® 

3.355*10"^ 

6ti 

\ 

10*"^ 

3. .14  *10"^ 

8.5<5*10"'^ 

7.76*10"® 

5.25*10"^ 

Another  approach  to  the  solution  for  a  magnetic  dipole 


(2,6)  and  (2,7)  gave  expressions  for  the  field  strengths  of 

magnetic  dipoles  whose  numerical  evaluation  is  somewhat 

difficult  and  can  only  be  simplified  by  restrictions  with 

respect  to  the  distance.  Besides,  in  the  expressions  of  (2,7), 

rather  complex  Kelvin  functions  occur  which  are  tabulated  only 

in  a  bounded  region  of  arguments.  In  order  to  avoid  this  we 

again  turn  to  the  integral  method  giving  an  expression  similar 

to  that  which  has  been  derived  in  a  thesis  of  the  Cornell 

University  [11],  In  the  latter  case,  however,  the  problem  of 

propagation  in  sea  water  is  dealt  with,  where  the  parameters  i 

(relative  dielectric  constant)  and  a  (conductivity)  assume  much 

higher  values ,  As  x  pe.  may  be  assumed  to  hold  for  the  quantity 

K  =  — ^  in  the  VLF  range,  the  approximation  in  the  above 
wto 

thesis  can  be  shown  to  be  vafid  also  in  this  case.  Furthermore, 


the  effect  of  the  ionosphere  shall  be  discussed  in  the  present 

r?  /  •  A 

2,  6,  Pig.  2,8,)  from  a  Hertzian  vector  n  =  (n  ,0,7i  ).: 
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*'  *9'  2 

H  -  k  71  +  grad  div  ti 


E  =  -iovi^i^rot  71 


(1) 


Let  the  time  dependence  be  neglected  and  let  it  be  exp(iwt) 
n  has  to  satisfy  the  differantial  equation 

2 


A  y~2^ 

^  71  +  k  71 


0 


U-^1 


The  condition  of  continuity  of  the  tangential  components  of 
H  and  E  make  the  following  relations  necessary  for  the 
"boundary  surface  2  =  0: 


dz 


dn 


Ax 


dz 


71 


Az’ 


‘  ~  !■', 


dis  n. 


E 


dis  7t, 


Ionosphere 


-h 


- 


Adr 


li;arth 


Pig.  2-5  a 


B  =’\(x^+ y^+  (z-h/  S'  .  Vx^+ y®  +  (z+h)^ 


If  a  homogeneous,  isotropic  and  infinitely  high  ionosphere  v/dth 
a  smooth  boundary  at  z  =  -h^  is  taken  into  accotint,  the 
equations  (;2)  are  to  he  supplemented  by  another  four  relations 
at  the  point  z  =  -h^  (the  sut scripts  E.  A,  J  refer  to  -,arth, 
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air,  and  ionosphere,  respectively). 


It  A  -  • 

Az  Jz 


A  Ax 


k^n- 
J  Jx 


div  n,  »  div  n_ 
A  J 


The  conductivity  of  air  was  put  zero  and  the  relative  dielectric 
constant  of  earth  was  neglected  as  compared  to  the  quantity 

H  ■  —  ,  We  thus  confine  ourselves  to  a  parameter  region  in 

which  the  displacement  current  may  he  neglected.  Thus,  the 
following  integrals  required  for  calculating  the  field  helow 


ground  are  obtained  by  means  of  the  well  known  integral 
representation  [S'J.  In  the  case  of  earth-air: 


®E  i  R 


E  A 


®  -Uj,(z+h)  ,  X 

.  V  JL  ^  e  ^  Ho^^  aVrUdX 

""ez  •  (l+-H)aj,coe(p  j  (n^+U  )(Ug-i>(U  ) 


In  the  case  of  earth-air-ionosphere: 


r  -ikgR  -tlCEjR' 
he  e 


^  R 

Im 


-%U+hO  A.y 

5  ®  Ho^^%ir^AdX 


-Up(z+nfA  /.A, 

,  ■  d  I  e  ^  Ho^^^(\r)XdX 


c 


Ruin,  2R.h^ 

E  A  A  1  . 

e  -  1 
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H 

4Tr 


are  valid. 


means  the  Hankel  function  of  the  first  kin'l  Re(UTn)y  0  is 

0  a 

assumed  to  guarantee  the  convergence  of  the  integrals  for 

(»a) 

z  -♦  OD  ,  and  Re  >0  is  assumed  so  that  the  factor  U^,  + 
in  the  denominator  of  does  not  vanish  in  infinity  of  the 
X-plane.  k  comparison  of  (4)  v'ith  (5)  shows  that  the  two 
representations  differ  by  the  factors  in  the  denominator  of 

2  2  1 

the  integrals  of  and  Since,  however,  =  h^(X  -k^)  ' 

is  real  on  the  path  of  integration  except  for  the  interval 
-k^<  X  <  +  k^,  and  is  a  large  number  (h~7*''0^  ^ti),  except  in 
the  somewhat  expanded  interval  A  -  k^,  v>'4-  k^  j  (v"  <  k^), 

we  may  put  tangh(U^h^ )  everywhere  equal  to  unity  by  approximation. 
Sirail-'^rly,.  we  put  F(X)  in  the  denominator  of  Tr^-,^  equal  to  zero. 
This  approximation  can  be  assumed  for  the  distance  region  in 

-ikjjR 

which  the  primary  excitation  — .  prevails,  since  the  other 

R 

terms  there  only  represent  corrections,  not  so  in  the  region 

z/r<^  1 ,  v;here  the  first  two  terms  in  cancel  out  and  the 

integral  expression  prevails.  For  an  exact  discussion,  another 

way  :has  to  be  chosen  for  this  region.  Since ,  however,  the 

distance  region  at  -Dresent  accesstbie  for  experimental  p  easure- 

3 

mentrs  below  ground  does  not  excedd  TG  m ,  we  certainly  may 
neglect  the  effect  of  the  ionosphere  there.  Hence,  we  apply 
expression  {4^  and  continue  calculating  the  ibtegrs-ls  U  and  V  j 


V 


7 

0  ' 


-Up(s+h) 


U,-  -  inU  , 


-00 


U 


. . .  1  inmiJ  ,  u<il  .lUVr-J  wliiiiir.  lUtUtfi'lUl 
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1  ^ 
2  in? 


°r  -Up(2+h)  (1),  , 

I  •  ^  8o  ^(Xr)XdX 

1  /w.wX/**  ^ 


V 

-c» 


E 


Substltutixig 


k  M  -  U, 

O  A 


f  - 


H  -  k  (z+h) 
0 


X  -  (7) 


we  obtain 


00 


-00 


r  -  ixM 


V  « 


JO  Jl^ 
2 


a> 


-00 


(N+M)(N-iwM) 


The  path  of  Integration  eomprisee  all  real  values  of  <};  ranging 
froa  -00  to  -foo.  In  the  following  calculations*  however*  it 
will  be  aore  convenient  to  shift  the  integration  path  slightly 
into  the  lower  semi-plane  at  the  point  (1  •  -1  (Fig*  2*b). 
Singular  points  occur  for 

I'  «  0 

M  -  0 

11  -  IkM  -0 


Thus*  branching  points  lie  at  <1  -  ±  f  (l  "i)  and  4*  •  ±  "*  • 


Since  the  latter  pair  of  points  would  lie  on  the  real  axis, 
it  is  more  convenient  to  write  down 


%  (1  -i;—) 


<Ci  -  t«l 


instead  of  putting  the  conductivity  oi  air  equal  to  zero  for 
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Pig.  2,3 


Thus ,  the  pole  lies  in  the  immediate  nei^'hborhood  of  Vg  • 
This  confirms  the  fact  (frequently  discussed  in  publications) 
that  the  contribution  of  the  pole  to  the  value  of  the  integral 
cannot  be  separated  from  the  total  expression.  Closing  the 
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path  of  integration  in  infinity  of  the  upper  aeni-plane,  the 

branching  sections  being  avoided,  we  include  Vg-and  in  the 

same  loop.  '  Since  for  large  ar^menta  behaves 

similarly  as  exp(i4>^)/(<i»^)^^^ »  the  contribution  of  the 

branching  section  may  be  ne,c^l*cted.  V/hen  integrating 

2  1  /2 

along  Cg,  we  must  consider  that  the  expression  K  *  (4;  +  in)  ' 
may  be  looked  upon  as  being  constant  owing  to  the  magnitude  of 
K  in  the  region  of  the  determina'tive  (j;-valuea.  Thus, 


-Hllx 


J  1  - 


1  - 


2\'n: 


is  obtained.  The  problem  is  thus  reduced  to  the  propagation 
between  two  points  of  the  earth's  surface,  z  »  0.  The  factor 
l/(l  -  ixM)  still  obtained  in  the  integral  is  thus  expanded: 


-  I’ix  M(i  -  -r-) 
ixM 


(1+7^  t-4r-  -K...) 

vix  M  i/ixM  )ix  M 


It  remains  to  be  shown  that  the  expansion  converges  for  all 
4rvalues  lying  outside  a  certain  circle  around  Vg.  The 
convergence  is  guaranteed  if 


j  'iiw  M  I 


j->»i  ,  i.e#  if  ^  . 


"»’e  show  that  this  is  the  Ci-.se  already : outside  a  circle  with 
the  radius  ^  ,  Vg  being  its  center: 


4;  »  -1+i|  +  ^  exp(i9) 


II  - 


“  ( i  -  ig)  +  \  (  ig  -  2  );  +  ^  I 


’Thus-j 


.2 


4  -  (1  -  ig) 


J  1  -  \I±H  Vl 

°2 


\  ig  -  2  + 


«v.  j 


00 


H 


\  V 
4L  n- 


(to) 


n=o 


,  1  .n+1  i  i4^'(>!-p)iM4- 

J  °  -,n;i 


‘.M 


is  valid.: 

In  the  case  n+1  =  2k(k=t,2,5,  •  •-• ) >  ^he  integrand  has  a.  pole  of 
the  order  (n  +  1  )/2  in  (Ij  =  -t.,  if  p:  is'  again  put  eqi'.al  to 
zero.  Henca, 


2vi+1 


(_!>_  .i2n  !  ^o 

IH  i  j  (^2  _  ^^p+1 


d=-t 


2ua  0  'S'' 


1  “  ix 


27ti  -1  ^  K' 

(dn)^  4 


can  he  v/riitten  down. 


In  the  case  of  integral  n,  v/e  a'*:'-  refer  to  a  paper  by  Somnei-- 
fcld  19}  according  te  v.'hlch 


f 


!  A 

t 

I  ^ 
1  y 

f* 


I 

.  I 


'  f 

4  ? 
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•!  2  2 


n 


-(4:  )"  ' 


I 


n-1  -i 


|t«  (n-1)!  j  Vo2  +o2 


00 


¥ 


£5 


1  w 


-1 


ViK5  5 


Thic  yields ; 

Ik  e  -  )/ix  H 

A  » 


X 


*■' "'■"’"(i;;?  f  (1  -i?)  - 


f<!  .  (i#  .ip 

- r  (1  +  1C  +  -  -  )  +  - — - pe  ^  + 

‘5  3  \|i« 


(11) 


IhhI’' 

,2  V^T  3 


for  the  first  terms  of  TJ. 

Confining  ourselves  to  ^-/x'^  1  v/hich  is  valid  in  a  \.'ide 
distance  region),  the  term  in  parenthesis  may  be  put  equal  to 
unity,  yielding 


U 


ik^exp(-i^,)  exp(-  \/iHH)  i  exp(-ikQr)  exp(- ^ix  k^(z+h)) 


? 


X 


X 


T'he  caicuiatioh  of  V  has  already  been  reduced  to  thai  of  U. 
HenCe> 


V  * 


ri  exp(-  l/ijk  ( z-fch^^ 


exp^^-ik-r) 

V  _ 


?<=l/r  ^  t 

.*r>V  ^ 


¥  iixi 


2x^  ix  k 


(12) 
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Ex  E 


exp(-ik^P\  exp.(  -ikgR  ’  )  iexp(-ik^r)  exp(- fix  k^Cz+h)) 


R 


R> 


(14) 


exp(“\^)k^(.z+h.)  exp(-k  r) 

’'Bz  -  ■=  — - j - ^d+i-k„r) 


The  first  term  in  denotes  the  direct  wave,  the  second  term 
the  wave  reflected  hy  the  earth's  surface.  The  third  term  as 
v.'ell  as  the  expression  for  can  he  interpreted  such  that 
part  of  the  energy  is  transferred  from  the  transmitter  through 
the  medium,  vertical  to  the  earth's  surface  along  v/hich  it 
propagates  above  ground  to  the  receiving  station  where  it  enters 
the  receiver.  Formally,  the  last  two  terms  correspond  to  the 
mode  denoted  as  secondary  wave  in  the  representation  of  (2,6)., 
since  the  exponent  of  e  in  IIj'  (2,6) 


i  k^(-nr  -  (z+h)'^/l+n^ ) 


may  also  he  written 


~lk  r  -  k  (s+h) 
o  '  o'  ' 


for  H 


This  representation  differs  from  that  in  (2,7)  the  folloT/ing 
points :  In  (2  , ,7)-,  the  integral  U  was  completely  neglected , 
this  is  admissible  only  if  transmitter  and  receiver  are  not 
too  clos:e  to  the  earth's  surface  and  a/r  (for  z  =  h)  is  not  too 
smalt  compared  with;  unity^,  since  in  both  cases  would  vanish 
alto  gather .  Furthermore ,=  the  integral  has  not  been  expressed 
by  the:  product  of  modified  Bessel  and  Hankel  functions,  but  by 
the  first  term  of  a  rapidi;;  convergent  series,  which  is  of  grf'st 
advantage  for  the  numerical  evaluation.  Thus,,  the  following 
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'expressions;  are  obtained  for  the  field  s.trength  components  t 

H,  :«p(.kgH):  51k/  ^ 

)  ^  2  ^^^E^  5x^  .\ 

-  H'  (■'e  -  HV  •  „.,.2  ^  „:,3  *  .p  ,  4 


_  exp(:-ik  r  )^  «  ik  1+kJx^ 

+  exp(-^  k^Xz+hO);  * >“7—  (k^  -  —  -  “2^ 

r 


Jik  x^  -  2 

+  i— 9,  ..  +  ) 

X  _  -A- 


exp(.ikj,R)-  .  k|  3ikg  ^ 

-V  ^  '"’Tv 


)  xy 


exp(;-ikgR»)  ^  i:  Jik^  5 

^  C  -  +  — r*  ~4  ) 

^  R'^  R’^  R'^ 


+  exp(-Yti<  kQ(iZH 

exp(-ik-  r  ) 

k^  3ik 

V  2  5 

x  r ' 

22  sxpX-ikgR) 

k|  5ikj, 

(-  + 
R  R^ 

“^r  )  X  (z-h> 
R^ 

,2  : 

exp^-il^R '  i 

Y 

+  •^^):  X  :(^z+) 
R*^ 

R'^  R=‘^ 

expj/ARO; 


ReHx  » 


+  R  sirt 


(Wf. 

r"**""  [ ^  cos it  +  31 4^in^ 


exp4^-?2i:C2H-h^  I:  „  . 

+,  ...I.; ^u--.^-;..-.-^.  i|s  ctoa  +  R  sin 
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Im  H 


exp'(*-AR^-i 


:R: 


B  cos  (^/J^R;)  «  A  sin  H 


exp|-i5^R:' ):  - 

■■'■  ~  '  t  ^  cos0iR  ')  -  C  sin  (^R'  )  J 


gxp(-A(zVh)) 


[  P  cos  (|)  -  E  sin  (|>  J 


In  the  main  direction-  of  emission  Hi  the  functions 

x|x=o 

A,-  B,  .  j-_  P  have  the  simple  form: 


B 


2  A^  - 
-  y 


R’ 


(17) 


4>  »  442+11)=  + 


In  order  to  determine  the  anguiar  dependence  of  the  value  of 
the  field  strength  in  the  plane  z  »  h,.  we  calculate.: 


at  r  =  R  =  const,  and  R'  =  const.  In  this  case,  the  following 
expressions  are  valid: 


H  =  A  +  Bx^  H  =-  Bxy  (I8) 

X  y 


exp4-lkj^^  - 

i  .  1,2  ^  n  -aJ 

Uxpl-lkjiR^l: 

;  ^  Jitkg 

i  * 

t  *  " 

R^ 

:1 

exp^(^-ikgR'4  : 

exp4^i 

kj,ii' )  i 

-4 

;  ^ 

;  R  -  H'- 

(4: 

■  R^'^ 

:R'^.  E‘^ 

I 

+  e 


xp$-$i 


^  *r»  - 


exp4-^^l i^jilc  fz+h ) -ik  r ' 


-I 

r  1 

' 

'  ,  2 

] 

i 

-  \ 
} 

»  ik  . 

^  ^ 

k4  5ik  i 

l^o  ■  r  ”  2^ 

T  ^ 

•*r 

=  ’  2  3  d 

r  T'^  r 

I 

1 
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•  1 2 

Thus,  for  (Hj;  we  ohtaln 

j  H);^  «  aI  +  |i2:  (iRe(j^^Re(B:)'  +  +  B^r  cos^^p 

or  [H#  -  +  E  .cos^<p  (1^> 

Re  *  i  t  real  part 

Im  imaginary  part,  A  Being  the  .conjugate  complex  value* 
iienoe,  j  in  a  polar  diagram  describes  the  curve  (19v  whose 
large  semi-axis  in  the  case  0’  lies  in  the  dipole  axis,  and 
in  the  case  B  <0  vertical  to  the  dipole  axis.  E>0  corresponds 
to  the  short  range  field.  For  , 

E  .  a|; 


T2r'^R 


9r 


10 


r®R.^" 


12R»‘'°  -  30r^R'^ 

“  8^,111 
r  R' 


>  0 


is  easily  found  for  E*  E<-G,  however,  is  valid  in  the  long- 
range  field.  Neglecting  the  terms  in  R  and  R'  (since  RA^R 
for  a  «  h  and  h/r  1 ):, 


- 

U  l^  -  :a|  ~  exp(-2d(z+h^  )  r^E 
r 


'  o®E 


exp(-A(z+h):) 


hence. 


-Er‘ 


and 


kVl-cos^(p=) 

®E  ^ - 2 - -  exp(  -2  A  (;z 


-r|- 

is  bbtained  for  z 


A  A 

»  JJK* 


Quatlitativelyv,  the  result  da  the  same"  as  in  dhe  /case  of  a  dtpole 
inaa  homogeneous;,  unbounded  medium:.  In  this  case 

r 

isF  valid  for  the  ahortTrange  f ieid,  whereas 


1h|:“  a  Gxp(r-2/ir)  *(l-cos^tp:) 

T 

is  valid  for  the  long-range  field. 

Thus,  a  quantitative  compa  ison  is  possihie.  The  ratio 
■between  the  emission  into  the  semi-space  that  into  the  full 
space  is 


for  the  long-range  field  (h/r  4  t).  Hence,  we  may  assume  that 
the  effect  of  the  surface  in  the  long-range  field  is  of 
considerable  importance.  This  is  also  confirmed  by  (2G),  since 
in  the  long-range  field  is  represented  only  by  the  term 
which  takes  the  existence  of  a  boundary  surface  into  account. 


z=h  =  200m 
(  =3kc 


J. 


i 
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Fig.  (2,3)  gives  the  course  of  the  field  strength  H  as 

xj  y— o 

dependent  on  the  distance  (up  to  r  =  5OO  m).  The  depth  of 

transmitter  and  receiver  is  600  m.  For  comparison,  the  field 

strength  of  a  dipole  in  a  homogeneous,  unbounded  space  was 

plotted  (dashed  curve).  Here,  an  effect  of  the  surface  is  not 

noted:  the  two  curves  overlap.  The  field  strength  decreases  as 

the  product increases  (x  are  measured  points).  Fig.  2,2 

and  Table  show  S’S  dependent  on  the  distance  in 

comparison  with  the  representation  according  to  (2,7)  and  the 

Hertzian  solution  for  the  unbounded  space  at  a  transmitter  ana 

receiver  depth  of  25O  m.  Fig.  2,2  shov;s  that  the  earth-air 

interface  hardly  has  any  effect  up  to  a  distance  of  1000  m. 

Table  2,6,  however,  shov/s  that  the  ground  wave  U  has  a  decisive 

4  5 

effect  already  from  r  =  10  onward,  at  r  =  10  it  yields  the 

10^^-fold  field  strength  in  the  case  o  =  mhos/ra,  and  the 

10^*^^ -fold  field  strength  in  the  case  o  =  5  *”10  ^  mhos/jn,  as 

compared  to  the  Hertzian  solution.  The  increase  in  field  strengt) 

5 

with  o  up  to  r  =  10  m  is  consistent  with  the  considerations  of 

['ll  (P‘  '15)*  The  dependence  on  a  is  reversed  as  the  distance 

increases,  the  field  strength  then  decreases  with  increasing  o. 

5 

The  propagation  in  the  range  r4  10  ra  is  strongest  according  to 
the  representation  in  (2,7)*  "^he  approximation  applied  in  (2,7) 
for  calculating  the  integral  expressions,  [in  (2,  )  the  follov/ing 
expressions  were  obtained; 


TJx 


exp(-ikgR)  exp(-ikgR‘) 
ag  (  R  HI  ) 
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=  a 
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exp(-ikpR’ ) 
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is,  however,  prohahly  not  suited  for  describing  the  propagation 
process  of  dipole  excitation,  for  the  following  reasons:  Owing 
to  the  boundary  conditions  (2)  at  z  -  0,  neglecting  the  integral 
U  in  successively  leads  to  the  conclusions: 


div 
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Introduction  of  the  expressions  for  7t^  of  (2,7)  into  the  last 
expression: 
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\z  “  I  exp(-Uj,^,)  J.,  (Xr)dXcos(p 


yields 


00 
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\z  ”  ®A  ®xp(U_^Z)  J^(Xr)d\cosq> 


yields 


az  az 


because  of  gp(X)  =  a  discrepancy  v/hich  cannot  be 


eliminated  by  the  approximation  used  for  because  of 


-"e  ^  \ 


or 


/  ]/x^ .  ki . 


At  di 


distances  of  r  ^  10^  ra,  R  and  r'  =  ^  r^  +  (z+h)^ = 


II 


for  z  =  h  ^  500  m. 


■become  approximately  independent  of  z 
Thus,  R  =  R'  =  r.  This  again  yields 
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k2 


?x?z 


[<>(oc)R„(p) 


If  the  asymptotic  represents-tions  for  the  modified  Bessel  or 
Hankel  functions  are  used 


\f2na 

K„(P)  =  e-P 

which  is  possible-,  since  for  ta|  =  A(il'  -  (a+h))  and 
|  p|  -  a(R'  +  (s+h)), 


I  a  j  ^  1  j  P  {  ^  1 
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holds  at  r>10  m  and  o>10  then  the  belov;  expressions  follov/ 
successively ; 
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exp(-ikjj(z+h)), 


ieS.  a  plane  v^ave  excited  in  z  =  -h  in  z-direction,  v;hose 

/  3 

amplitude  decreases  as  1/r  .  An  exact  valuation  of  the  theory 
described  in  (2,7)  v/ill  be  given  later. 
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III.  Comparison  of  theory  and  experiment 

The  measurements  described  in  part  I  v;ere  evaluated  from  the 
viewpoints  of  two  different  theories,  one  being  the  theory  for 
a  dipole  in  an  unbounded  homogeneous  medium,  and  the  second 
being  an  elaboi’ate  theory  developed  in  chapter  2  and  also  in 
Unfortunately  it  has  not  been  possible  to  evaluate  this 
theory  accurately  and  in  detail,  as  for  this  purpose  an 
electronic  computer  would  be  necessary.  We  shall  make  these 
calculations  in  the  near  future  ,  as  a  computer  is  going  to 
be  set  up  in  Innsbruck  very  soon.  So  far,  we  could  compute 
only  some  special  measurements  with  considerable  simplifica¬ 
tions. 

The  measurement  in  the  German  mine  Fusseberg,  dealt  vvlth  in 
chapter  1,7,  seemed  to  be  well  suited  for  studying  the  effect 
of  the  earth's  surface  on  wave  propagation  over  large  distances, 
because  the  geometrical  conditions  are  favorable  there. 

The  points  of  measurement,  the  curves  for  a  pure  dipole,  and 
those  according  to  V'j  and  chapter  2,8  were  plotted  in 
Pig.  2,5,  but  a  statemeni  on  the  exact  course  of  the  experimental 
curve  cannot  be  made  yet  because  the  distance  is  too  small 
(only  500  m) .  In  this  distance  range,  at  a  conductivity 
0  =  T0”^52  '^m~^  which  in  this  mine  can  be  assumed  to  be 
mean  value,  the  two  above  theoretical  curves  are  practically 
identical . 

The  conditiens  of  our  measurements  in  Schwaz  are  similar  -  here, 
the  measurement  reaches  1200  m,  but  the  effect  of  the  earth's 
surface  cannot  be  proved  witr  certainty  because  of  the  low 


conductivity  and  the  great  depth  (Fig.  5»2), 

For  our  measurements  in  Schwaz,  the  surface  can  no  longer  he 
assumed  as  plane  (as,  e.g.  in  Tlfsseberg). 

Table  5»1  gives  the  values  calculated  according  to  the  theory 
of  a  dipole  and  according  to  f?!* 


The  transmitter  described  in  chapter  1 which  was  used  for 
measuring  over  distances  up  to  3 *500  m  in  the  St.  Gertraudi 
mine  showed  good  results  (chapter  i»9)» 

The  geometrical  conditions  of  these  measurements  are  given  in 
Fig.  3  »5  and  the  idealization  necessary  for  the  evaluation 
is  shown  in  Fig.  5 >4 


T— .r-  l  V 


^3 
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-  ^  ’  V 

""  '  earth's  surface 


EW  profile 


(x)  indicates  the  position  of  the  transmitter  and  (:o)  indicates 
the  individuai  points  measured  below  ground  as  well  as  on  the 


earth ’ s  surface . 
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Fig.  3»4  shows  that  the  depth  whi.ch  had  to  he  introduced  into 
the  calculation  decreases  fbr  the  individual  points  of 
measurement,  Despite  the  considerable  idealization  of  the 
present  problem,  the  course  given  in  Table  3»2  and  Fig.  3>1 
shows  that  the  points  of  measurement  at  distances  of  more 
than  1.5  km  are  much  better  described  by  the  new  and  improved 
theory  (allowing  for  the  sky  wave  and  gi’ound  wave)  than  by 
the  original  dipole  theory. 


It  will  be  the  aim  of  further  measurements  with  a  more 
pov/erful  transmitter  antenna  which  is  now  being  built  to  study 
the  propagation  over  very  great  distances  in  detail. 

In  this  case ,  only  the  ground  wavs  yields  a  considerable 
contribution  to  the  propagation  mechanism. 


Table  3,1 


d  =  5.10"'^  >  transmitter  -  earth’s  surface  1000  mj  f  =  3  kc. 


Dipole  1^1  j 

tn6  [tJ  j 

1  measurement  |^Wb/m‘ 

200 

-1 0 

2480*10 

2670.10-IC 

1 

3,8.10*'''2 

1 

1 

516 

626 

657 

j 

405 

306 

514 

'  1 ,5 

620 

81 

92 

0,38 

740  ;| 

45 

51 

0,21 

850  j 

51 

34,2 

0,13 

910  I 

23,8 

27,8 

0,095 

1000 

17 

20,5 

0,07 

1024 

16,6 

18,7 

0,07 

11 80 

10,3 

10,7 

0,042 

1275 

7,87 

6  ,  08 

0,C55 

Ill -4 


Table  5,2 

St.  Gertraudi;  h  =  400  m,  a  =  ra"\  f  =  5  kc,  6  =  50^ 


r  !_.]  I 

1  Ih{[7]-10'^^ 

measurement  jjivJ 

1 

1000 

11,5 

0,45 

1100 

8,81 

0,445 

1200 

6,85 

0,55 

1400 

4,2^  ; 

0,262 

1550 

5,55 

0,548 

1800  1 

! 

1 

2,20 

0,151 

2000 

1,5 

0,088 

2200 

1 

1  ,1 2 

0,197 

5000 

0,507 

0,058 

5500  ; 

0,408 

0,01 97 

1 

4000  j 

I  0,150 

1 
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The  publications  (I)  through  (4)  as  well  as  (14)  pr^'ved  to  be  of 
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principal  importanci .  The  methods  discussed  therein  were  used 
in  Technical  Notes  Nos.  4  ’to  6.  An  exact  study  of  (7)  made  us 
acquainted  with  fundamental  methods  for  the  treatment  of 
transmission  problems  in  general,  and  especially  those  in  the 
VLF  range.  The  Mode  Theory  oxolained  in  (5),  (6),  (7),  (9), 
and  (l0)  as  well  as  the  methods  of  VLF  transmission  in  the 
atmosphere  between  the  earth's  surface  and  the  ionosphere 
explained  in  (l5)-»  and  (19)  through  (24),  where  the 

w«,,v~  are  considered  to  be  excited  in  the  atmosphere,  v/ill  be 
important  for  further  work  because  the  theory  may  be  extended  to 
excitation  of  the  modes  beyond  the  atmosphere;  interesting 
results  are  expected  to  be  gained  by  this  extension,  A 
continuation  of  the  above  studies  is  planned  but  is  bound  to 
the  use  of  an  electronic  compucer  which  will  be  available  at 
the  beginning  of  1964.  As  yet,  publications  (17)  and  (18) 
seem  to  be  of  little  interest  for  our  work.  Publications  (25) 
through  (30)  as  well  as  (8)  deal  with  antenna  problems. 

In  this  connection  it  will  be  necessary  to  study  an  even  larger 
number  of  publications. 

The  comprehensive  report  (62)  which  we  obtained  o'Jst  before  going 
to  press  naturally  was  of  special  interest.  It  showed  cn  the 
one  hand,  that  the  way  we  chose  independently  for  describing 
the  propag^ation  of  VLF  waves,  has  also  been  used  by  other 
scientists-.  It  proves,  on  the  other  hand,  that  a  nuntber  of 
questions  must  be  answered  before  the  problem  can  be  looked 
upon  as  being  solved,  although  some  problems  can  be  solved  by 


mcdify;ing  the  existing  theo^i^y  of  i?he  propagation  and 
excitatrion  ol  H,F  -wamar  in^  the  atmosphere  hetween  earth  and 
ionosphere.  This  holds  both  in  theory  and  practice, 
especially  for  the  antenna  problem.  A  considerable  improvement 
of  the  efficiency  of  VLF  antennas  would  be  a  great 
advance  in  this  field,  which  has  not  yet  been  made  despite 

a 

great  efforts.  For  further  studies  of  the  excitation 
problem:,  emphasis  must  be  laid  on  improving  the  efficiency 
of  VLF  antennas. 

(32)  solves  the  problem  of  the  propagation  of  electromagnetic 
waves  in  the  interior  of  the  earth  in  three  different  ways, 
depending  on  the  geological  structure  and  on  the  values  of 
the  electrical  parameter  of  the  studied  region: 

(:a)  by  the  so-called  UP-OVER  Allf  DOWN  modes  at  a  small 
depth  of  transmitter  and  receiver,  and  at  high  conductivity 
losses,  a  form  of  propagation  R.  E.  Moore  described  in  this 
thesis. 

Here,  the  direct  wave  and  the  wave  reflected  from  the  earth's 
surface  are  neglected  because  of  high  exponential  damping. 
Hence,  this  theory  can  describe  propagation  at  high 
conductivities  (-approximaytely  that  of  sea  water)  . 

(?b^  Dy  the  so-cailed  SHORT  DISTANCE  MODES  at  medium  ohmic 
iosses;.  The  field  p^t  the  point  of  reception  is  composed  of 
a  direct  wave ,,  a  refiented  v/ave  ,  and  a  ground  ware  (  we  have 
already  deacribed  ibla  f<?rin- of- propagation  in  j  7^ )  • 


-  By  th  guide  theory  ;a-t  -smaOO:  ohmic  la^sestt  as  i-s 

the  case  in  propagation:  In  the  ^atmosphere  -betv<een=  earth,  and 
ionosphere;.  The- geoiogical  structure  of  the  earth  was 


assumed  to  consists  of  layers  of  difierent  conductivity  (^st  layer 
oi  comparatively  poor  conductivity  between  two  iayers  of  good 


c  ondu  cti  V  iiy ): . 

(br)  has  already  been  worked  out  in;  ^7^  and  has  been  adepted  in 
the  present  report,  (a)  and  (c;)  are  of  less  interest  for  our 
work,  because  the  geological  conditions  are  not  appropriate. 

Our  measuremente,  are  mainly  carried:  out  with  the  conditions 


shown  in  Pig.=  i,/-*  Therefore  we  are  intended  to  set  up  a  theory 
of  propagation  based  on  this  st  ructure ,  which,  will  better 


correspond  to  the  actual  conditions:. 


Ionosphere 


atmosphere 


ground -waier  ^  -.'jf-Jm 
„  level  _ J _ 


Another  problem  will  be  to  study  the  propagation  in  a  layer 
with  a  conductivity  <j  =  ^  mhos/m  making  allowances  for  the 

refiactioir  from;  the  aground  and;  hy  the  grounds- water  levefl.^, 
as  well  as  the  influence  of  the  i*n# sphere  which  will  exist 
at  least  for;  large  dietanoes  ^nurre  than  5®0Q  m^l^tth  thh^^^ 
transmitter  and:  rejceiver  at  comparatively  small  depths 


The  present  paper  deals  with  measurements  concerning  the  propagc-- 
tion  of  VLF  waves  through  inhomogeneous  rock,  the  influence  of 
the  earth’s  surface  being  taken  into  account.  These  measurements 
are  compared  with  results  of  theoretical  works.  Furthermore, 
conductivity  measurements  of  rock  in  mines  are  discussed  and 
compared  with  the  results  gaiued  from  propagaticn  measurements. 
The  results  have  been  found  to  agree  well.- 
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